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Abstract 

The hydrodynamics of Newtonian fluids has been the subject of a tremendous amount of work over the 
past eighty years, both in physics and mathematics. Sadly, however, a mutual feeling of incomprehension 
has often hindered scientific contacts. 

This article provides a dictionary that allows mathematicians (including the author) to define and 
study the spectral properties of Kolmogorov-Obukov turbulence in a simple deterministic manner. In 
other words, this approach fits turbulence into the mathematical framework of studying the qualitative 
properties of solutions of PDEs, independently from any a-priori model of the structure of the flow. 

To check that this approach is correct, this article proves some of the classical statements that can be 
found in physics textbooks. This is followed by an investigation of the compatibility between turbulence 
and the smoothness of solutions of Navier-Stokes in 3D, which was the initial motivation of this study. 

The simplest model of a Newtonian fluid is an incompressible flow evolving freely with constant density 
and temperature. Let us therefore consider the incompressible Navier-Stokes system on M + x Q with either 
HCE 3 orO = T 3 : 

' d t u - vAu + (u ■ V)u = -VII 

< divu = (1) 
u\ t=0 = n , u\ dn = 0. 

Here u represents the velocity field and n = p/p, where p is the pressure and p is the density assumed to be 
constant and v > is the kinematic viscosity. This equation has weak solutions (called Leray solutions [50]) 
in the Leray space 

u G £(fi) = L oo (R + ;L 2 (^))nL 2 (R + ;ij 1 (0)), n G L^(R + x fi). 

Pressure can be computed by solving —AIT = Tr(*Vii • Vu). In general, the smoothness of (u, II) is an 
open problem that has been extensively studied. Among historical landmarks regarding smoothness, one 
must cite the works [60], [32] and [39] for the point of view of partial differential equations and [14] for 
an approach based on geometric-measure theory. To get a more comprehensive survey of what is currently 
known about (1), one should check e.g. [48], [54], [65] and the references therein. 

Since the seminal works of A.N. Kolmogorov [40], [41], [42], [43] and A.M. Obukhoff [57], a vast amount 
of effort has been put into understanding turbulence. In physics, one should definitely quote [2] and [31] as 
major reference handbooks. Personally, I was very impressed by experiment [53], in which immersed floats 
equipped with GPS devices were allowed to drift in a Canadian river; the speed could be measured directly 
and a sufficient amount of data could be gathered to check the spectral 5/3 law (see Prop. 5 below) with 
striking precision. I discovered the point of view of engineers in [25] and was pleased to realized that they 
pay great attention to mathematical rigor because disregarding the divergence of an integral can trigger a 
catastrophe in real life. The following books and articles helped me acquire the experimental background 
necessary to write this article [8], [6], [7], [10], [56], [63], [35], [37]. In mathematics, the question of turbulence 



often raises cynical reactions. However, the books [30], [52] and the following works were a valuable source 
of inspiration : [22], [46], [23], [45], [1], [29], [24], [20] (by publication date). 

On the question of whether the spectral properties of turbulence are compatible with Navier-Stokes, 
there is a very interesting recent paper [5] dealing with weak solutions. The authors obtain compatibility 
conditions that they themselves qualify as being reasonably satisfied (the upper-bounds on the inertial 
range largely exceeds the range predicted by physics and observed experimentally). The present article is an 
independent work, though my motivations are similar. Starting from a more precise definition of turbulence 
makes it possible to recover the real inertial range. Later on I focus on smooth solutions and prove the 
exponential decay of the spectrum, which will lead to much stronger restrictions that still allow smooth 
turbulence to exist. 

Structure of the article, ideas and main results. 

Section §1 contains useful definitions and notations. In particular §1.6 provides a deterministic definition 
of K41-turbulence and a dictionary to translate physics claims into mathematical statements. Section §2 
checks that it indeed leads to the classical properties of the inertial range and of the energy spectrum that 
one can find in physics textbooks. 

The next step is to investigate a-priori bounds of the energy spectrum respectively for low frequencies 
in §3 and high frequencies in §4. The low-frequency bounds happens to answer a physical conjecture on what 
triggers Batchelor's and Saffman's spectra ; the answer is strongly connected with the spatial localization 
and temporal decay properties that have been extensively studied by mathematicians in the past decade. 

Section §5 contains the discussion of whether solutions of (1) in 3D can be turbulent and what conditions 
must be satisfied. A general restriction applies to the time range on which averages are taken. For smooth 
solutions, one can prove a lower bound of the fluctuation between the dissipation rate and its average (a 
physical phenomenon known as intermittency) . There is also a new formula relating the analyticity radius 
to the size of the finest scales in the inertial range. 

As some other important statements of physics textbooks have not yet been rigorously established, 
section §6 collects some open problems and hints on how they could be tackled. 

The first main idea carried by this paper is that qualitative properties of turbulent flows can be studied 
with the deterministic tools of PDEs. A probabilistic approach might still be necessary later on to prove 
that "most" flows are turbulent, but that problem should be addressed separately. 

The second idea is that turbulence is not based on the failure of smoothness because one can prove it to 
be compatible with the analytic regularity of solutions. Turbulence is a specific mathematical structure of 
solutions that local regularity methods fail to capture even though they might be troublingly close (see §6.8). 
In all likelihood, turbulence will prove to be the key to understanding global smoothness. 

The third and more concrete contribution of this article is the introduction of the "volume" function 
Vol (it; [To,Ti]) at the foundations of the theory. This quantity describes the large scales involved in turbu- 
lence. Later on, the new time scale T(uq',uj) is shown to be characteristic of time intervals on which free 
turbulence can be observed. Their ubiquitous nature in this article suggests that they should be given some 
attention. 

Finally, this article sheds light on some subtleties related to the definition of the energy spectrum in 
the discrete case. Luckily enough, spectral computations could be carried out explicitly on T 3 but the 
generalization to other domains should be done very carefully as it might be responsible for a substantial 
part of the troubles of "real life" turbulence. 
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1 Definitions and basic properties 

Let us recall some mathematical notations and physical definitions regarding (1). 

1.1 Kinetic energy and energy spectrum 
The kinetic energy at time t > is 

E(t) = p\\u(t)\\ 2 L2 . (2) 

The energy spectrum represents the contribution to the total kinetic energy of the frequency K. It is 
defined rigorously by the spectral resolution of the Stokes operator A = —PA which is a positive self-adjoint 
operator with domain D(A) = {«£ H 2 (Q) n Hq(Q) ; divii = 0} if fi is smooth (for non-smooth domains, 
see e.g. [61, p. 7, p. 128] or [33]). Here P denotes the orthogonal projection on divergence-free vector fields 
in L 2 (Q). 
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1.1.1 Case of a continuous Stokes spectrum 

On f2 = R 3 one has a(A 1 / 2 ) = [0, oo). Thus the (isotropic) energy spectrum of a function u G £(M 3 ) is 
defined by : 



E*(K, t) 



d 
dK 



X 



'a 1 ' 2 ' 

K 



u(t, 



with 



L 2 , 



r xeC°°(I 
X(r) = 1 



r < 1/2, (3a) 
r > 2. 



We use a smooth cut-off to ensure that E* exists for any u(t) £ L 2 (M 3 ). The non smooth cut-off is defined 
as the limit in the distribution sense x 2 ( r ) - 1 lr<i- The spectrum satisfies the fundamental property : 



/•oo 

£(£) = / E* (K, t) dK 
Jo 



(3b) 



The projector P commutes with derivations, which allows one to compute the spectrum explicitly in Fourier 
variables (see e.g. [19, pp. 38-40]). 

Proposition 1 Let us compute the Fourier transform with u(t,£) = f R3 e~ tx '^u(t,x)dx. 
1. If u is a divergence-free vector field in R 3 : 



K J , ~ v ''*" (2vr) 3 

where ip is a positive smooth bump function, supported on [2 _1 , 2] and such that J °° ip( r )^r = 1- 
2. In the limit of non smooth cut-off, one has ip(r) — 5 r= \ and 



(4) 



E*(K,t) ->Efl(K,t) = (2tt)- 3 P K 2 [ \u(t,K0)f 

dcf J S 2 



d-d. 



(5) 



Conversely, if one defines the "experimental" value of E^(K, t) as the average over a spherical shell of 
relative amplitude 5 e]0, 1[, one finds 

r(l+8)K 



- / \E^K,t)\ 2 dn = E*(K,t) 

2dK J (1 -s)K 



'{1-5)K 



for ip(r) = ^l[i_5,i +( 5](r) and ip(r) -± S r=1 as S ->• 0. 

Proof. Applying the spectral theorem, one has A = J °° AdP\ with P\ = F o P^ where P^ is the spectral 
projector associated to the Dirichlet operator (—A)d on Q. For M 3 and T 3 , one has Pf = J^oPfoJ 
where is the Fourier transform and respectively (£ £ R 3 and fcsZ 3 ) : 



or P\ = l|fep<A- 



In both cases, one has P o PP = PP o P hence 



X 



A 1 / 2 
K 



u = F- 1 o 



X 



\ 1 ' 2 
K 



dP x oT(u) =P~ l 



X 



K 



Fu 



If divu = 0, one has Pu = u. Then (4) follows from Parseval identity 1 1 1 1 1, 2 (IR 3 ) = i^n) 3 ||w||^2. The rest of 
the statement is obvious with ijj(r) = —2rx{r)x'(r). ■ 
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Remarks 

1. To unify notations with the case of a discrete Stokes spectrum, let us state (3b) as 



E(t) 



E\K, t) dfi(K) 



(6) 



where is defined by (5) and jx is the Lebesgue measure on a{A 1 ^ 2 ) = [0, oo). Note that by Fubini's 
Theorem, the spectrum E^(K,t) exists in L^IL^) for any u E L 2 (R 3 ). 

2. In the following, one should not rely on any other norm than ||^||^i < 2 because other L p norms of ip 
are unbounded in the limit ifj(r) — 6 r =i. Estimates like E*(K,t) < \\ip\\ Lx should be disregarded 
as empty of any physical meaning and because they do not correspond to any property of . 

1.1.2 Case of a discrete Stokes spectrum 

On T 3 = R 3 / (LZ) 3 the spectrum of the Sokes operator is discrete cr(A) = {K 2 ; K E E} where 

E = a(A 1 ^ 2 ) = {|k| ; k E (2nL~ 1 Zf} = {27rL _1 Vn; n E D 3 } (7) 

and D3 is the set of integers that are the sum of three squares. Let us denote by E* = E\{0}. By the 
Gauss-Legendre three-squares theorem : 

□ 3 = {nEN, n/4 p (8g + 7)}. (8) 

Writing an analogue to formula (3a) must be done carefully in order to cope with the discrete differentiation. 

Let us denote by u(t, k) = f T3 e~ lk ' x u(t, x)dx the k th Fourier coefficient of u(t) and by L = Vo^T 3 ) 1 / 3 
the characteristic length of T 3 . 

Definition On T 3 the energy spectrum is defined by : 

VKeE*, Efl(K,t) = (2ir)- 2 p(£\ £ |n(t,k)| 2 . 



(9) 



ke(27rL- 1 7^ 3 
\k\=K 



One defines the following measure on E* 



For 5 <C K, the corresponding average value on the spherical shell T,$(K) = |k E E* ; ^ 



< 



K 



E* s (K,t) = " MK) 



J-Es(K) 



2.l- y: 



2tt 
kL 



(10) 
}*: 

(11) 



\kL J 



One will call E* s {K,t) the "experimental" value ofE^(K,t). 

Proposition 2 The following statements hold. 

1. To compute the total energy, formula (3b) is replaced by : 



E(t)=(f pu (x)dx] +[ E\tt,t)dK 

\iT3 / J 1 jM1/2)\{0} 



(12) 



f pu (x)dx) + 2-KL' 1 V ( — )E\K,t). 
1 ■ ' Kex« 



( 2vr 
\KL 
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2. There exist /3 e]0, 1[ and C > 1 such that 

0<S<(3K => C- l El{K,t)<^-\^ Y. \Ht,k)\ 2 \ <CE* 5 (K,t) (13) 




where S$(K) = j& € (27rL 1 Z) 3 ; I ^ — 1 < ^| is the spherical shell of frequencies K ±5. 

Identity (13) is crucial to match the theory to real-world experiments. Indeed, the sum of the squares of 
Fourier coefficients on spherical shells is the energy spectrum of all numerical and physical experiments. 
Therefore, the universal behavior observed for the Eg of turbulent flows can be addressed mathematically 
by investigating the corresponding property on E^ . 

Remarks 

1. One cannot emphasize enough that (13) does not hold for any other normalization than (9)-(12). For 
example, using the analogy with (5) one could be tempted to replace (9) by : 

(2vr)-V^ 2 E K^)! 2 

Lk£2nZ 3 
\k\=K 

This choice would lead to a catastrophe. First, the averages on spherical shells would be equivalent to 

I 772 E 




which is not the usual normalization of experimental spectra. If this fact remains unnoticed and 
one develops the rest of the theory, then in Theorem 6 one would not recover the usual Kolmogorov 
dissipation frequency (defined below) but instead K' d = - . Physics textbooks usually dodge 

the subtlety of (9)-(12) ; at best, they use (5) to define E*(K,t) in the case of M 3 and rely on (13) for 
any practical purposes. This subtlety is however of great practical importance as most experimental 
data is obtained in a situation where the spectrum is discrete. . . 

2. The key to this computation is the asymptotic of S* i.e. of eigenvalues counted without multiplicity. 
Let us reorder the Stokes spectrum in an increasing sequence a(A) = {K 2 ; j G N} with Kj < Kj + \. 
One can prove the existence of C, C > 1 such that : 

^ < K] +l - K] < g or equivalents - £ * < K j+1 - Kj < (14) 

This asymptotic is responsible for the fact that E^ (or more generally any function on o{A 1 / 2 )) has 
a different normalization than its shell averages Eg. Describing which domains satisfy (14) on the 
pertinent range of frequencies is a widely open problem whose answer might actually provide a hint 
as to why physicists have trouble unifying the description of all turbulences. Conversely, it might 
also explain why some common patterns have been found for each of the "families" of turbulence 
{e.g. grid turbulence, wake turbulence, jet turbulence, boundary layer turbulence or turbulence in 
bounded domains). 

3. The experimental spectrum is defined by averages on spherical shells of frequencies K ±5 with 5 <C K. 
For dyadic shells i.e. frequencies (1 ± e)K the corresponding average is 




E: K (K,t)^^[j- 3 e !«(*.■ 



-e<\k\/K<l+e 
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Proof. Parseval indentity 1 1 ^ 1 1 ^2 (ir 3 ) = ^ 3 2~^l^(^)| 2 dictates the compatibility between formulas (9) and 
(12). Note that (1) implies / pu(t,x)dx = I pu^{x)dx but the corresponding energy is not seen by the 
spectrum. Parseval identity also ensures that the numerator of E$(K, t) is 

To prove (13) one has only to check that the denominator 



/2tt 



2vr 



1 



n'en 3 



where 



K 



2vr 

T 



n 



is equivalent to 5, which is an easy exercise in number theory (see the proof of Theorem 6 below, where 
a similar computation is fully detailed). The numerical value of the constant is illustrated in Figure 1. 
The key ingredient is that N\d3 contains only integers n = 0, 4 or 7 mod 8, which in turn ensures that 
Xme[i N]nn 3 n ~ S an d En=i n ~ s are equivalent up to a numerical factor. This exact same property of con- 
gruence modulo 8 also implies that 



n G D 3 =^ □ 3 n{n + l,n + 2,n + 3}^0. 
The asymptotic (14) of the Stokes spectrum follows immediately as y/n~+j — v'" 



(15) 



-J- + O 




8000 10000 



Figure 1: 



Plot of ^ V 4=7 for n' G D 3 such that 

1 



n < — \/n. 
- 10 



On this range, one has 



10 X 1 

3.3 < — 2^— J (>.2 



The numerical value of the asymptotic equivalent is 
3.335. The corresponding sum without the restriction 
n' G O3 is equivalent to 4 as n-> 00. 



1.2 Dissipation of kinetic energy 

The dissipation rate at time t is defined by : 

s(t) = 2pu\\Vu(t)\\ 2 L2 . 
One can check immediately in Fourier variables that : 

-e(t) <2v \ K 2 E*(K, t) dK < 2e(t) on 



- Jo 

e(t) = 2v I K 2 E\K,t)diJL(K) 



on T d 



(16) 

(17a) 
(17b) 



Note that 2 \\Vu(t)\\ 2 L2 



L 2 if u is a square integrable divergence-free vector field on M. 3 or T 3 with 



vorticity u = (diUj — djUi)i<ij<3. Because of the so called "stretching term" in the right-hand side of 

dtoj — vAco + {u ■ V)w = (u) ■ V)u 
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no good a-priori estimate of uj or e is known (check [15] for partial results). One has only : 

e(t) = e(0) - 2pu 2 f ||Vw(r)||^ dr - ApuY] f I (d iUj )(d kUi )(d kUj ). (18) 
Jo lJ k Jo Jn 

Gronwal inequality provides exponential estimates like : e{t) < Ce(0)exp ^ J* || Vu(t')\\ Loo dt'^j. 



(19a) 



Global balance of energy. If u is a smooth solution of (1), one has 

dE 

e ~ ~~dT 

which justifies the name of "energy dissipation rate" for e. As e depends only on co it means that the 
dissipation of kinetic energy occurs exclusively through vortex structures. For Leray solutions, one has only 
for a.e. t > t' : 

E(t) < E(t') - [ e(r)dT. (19b) 
J? 

As observed in [26] the possible lack of smoothness (i.e. a strict inequality in (19b)) would mean that an 
extraordinary dissipation has occured between t and t! . The discussion of whether (1) would remain a good 
physical model on such a [t, t'] (or what model should replace it) is beyond the scope of this article. 

1.3 Scaling transformation to work per unit of mass 

One can construct two families of transformations that preserve (1) without changing the kinematic 
viscosity v. 

• If (u, IT) is a solution of (1), then so is : 

VA > 0, u x {t, x) = Xu(X 2 t, Xx), U x (t, x) = X 2 U(X 2 t, Xx). (20a) 

In the case of T 3 , the new domain becomes T x = M 3 /(A _1 Z) 3 and to ensure the conservation of the 
total mass, the new density must be : 

px = X 3 p. (20b) 

• (u,p,p) I—?- (u, [ip, pip) for p > is another family of solutions. This transformation means that the 
mass of each particle is multiplied by a factor p without changing the number of particles in the fluid 1 . 
In the case of T 3 , this transformation would obviously change the total mass by a factor p. 

Definition In the following, one shall work "per unit of mass" i.e. for a given number of molecules. For 
a bounded domain Q or T 3 , it means that one choses p = I / pdx I and is left with 



p = Vol(n)- 1 . (21a) 

For M 3 , one will extend this by convention by letting : 

p = 1 [length]' 3 . (21b) 

Once one works per unit of mass, one cannot apply the second family of transformations anymore. The first 
one (20a)-(20b) remains admissible with a dimensionless parameter A > 0. It acts in the following way on 
the energy-related quantities : 

E x (t) = X 2 E(X 2 t), E{(K,t) = XE*(j,X 2 t ] j, e x (t) = X 4 e(X 2 t) (22) 

and ||u A || y = \\u\\ Y for Y = L°°(R + ; X) and e.g. X = L 3 (R 3 ), H 1 / 2 ^ 3 ) or BMO^iM. 3 ) or any of the so 
called "scaling invariant" function spaces. 

Like replacing a hydrogen flow by a helium flow with the same velocity field and assuming that all other physical properties, 
including v, will remain identical. However these particular gases are compressible so (1) does not describe them properly. 
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Figure 2: Real-world illustration - kinematic viscosity of 
liquid water as a function of temperature (see e.g. [44] 
and the references therein) ; the order of magnitude is 



v ~ 1(T 6 m 2 s -1 . 



Physical dimensions. Let's recall that v = [length] 2 - [time] -1 is the kinematic viscosity. For liquid water, 
the value of the kinematic viscosity is illustrated by Figure 2. One can easily check that the natural scaling 
of energy per unit of mass is 

E(t) = [length] 2 • [time] -2 . 

It follows that E*(t,K) = [length] 3 • [time] -2 and e(t) = [length] 2 • [time] -3 . With the notations of Proposi- 
tion 1, the Fourier transform is u(t, •) = [length] 4 • [time] -1 on both M 3 and T 3 , the modulus of frequencies 
is K = [length] -1 and the density is p = [length] -3 . From a mathematical point of view, working with 
physical dimensions is equivalent to checking that each identity is scale-invariant under (20). 



1.4 Time averages and intermittency 

The function E*(K,t) and e(t) are to some extent accessible to the experiments of fluid dynamics. And 
it is a common observation that even if those functions fluctuate a lot, time averages (over proper time 
intervals) display universal behaviors. This observation is the experimental essence of turbulence theory. 
Time averages of a function f(t, 9) on [To, T\] x will be defined by : 

f( ) = X r /(*> °) dt with A = Ti - T . (23) 

We will also denote / by (/) if the expression of / is too large and makes the first notation ambiguous. 

Definition Once the time interval [To,Tl] is given, one defines the mean energy E, the mean energy 
spectrum E*(K) and the mean energy dissipation rate e according to (23). 



Substantial fluctuations of e{t) away from e are an interesting phenomenon known as intermittency. 
More subtle definitions are possible ; one could call this one "temporal intermittency" to differentiate it 
from "spatial intermittency" that deals with substantial spatial fluctuations, either at large or small scales. 
The following result provides a simple way to detect intermittency by comparing the average energy E to 
the linear interpolation between the initial and final energy. We will need this statement for Theorem 16, 
which establishes a subtle relationship between intermittency and the smoothness of turbulent flows. 



Proposition 3Ifu is a smooth solution of (1) on [To,Tl] then : 

Eq — Ei 



A 



(24) 



with Ei = E(Ti), A = T 1 -T and 



E 



E + E x 



< I \e(t) - e\dt. 

'T 



(25) 
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Proof. As u is smooth, the balance of energy reads E(t) = E(t') — J*, e(r)dr hence e = E ^ T °^ A E( ~ Tl ^ ■ Let 
us now integrate this relation on [To, To + t] for t G [0, A] : 

E = 1 / jB(T ° + ^ dt = E ^ ~ 1 j T 1 (T ° + A " r)e(r) dtdT 

hence : 

- _ Eo^E, = -_ E{To) + A_e = 1_ ^ (Ti _ r)(£(r) _ _ } 
One concludes using the L°° * L 1 — > L°° convolution property. ■ 

1.5 Average "volume" of a function 

In naive terms, one can describe Vol (it; [To, TJ) as an intrinsic measure of {x ; \u(t, x)\ > e} for "adequate" 
e and proper time average, i.e. the average volume of the region where u is most intense. 

Definition For any measurable function u(t,x) G L 2 ([To,Ti] x fi), one defines the average volume 
occupied by u on [To,Tl] by : 

/ll ll ^ \ 

VaL(u;[r ,r 1 ])= , f (n) € 1+ (26) 
where (•) refers to time averages defined by (23). 

The following sections will show the relevance of (26) for turbulence. Definition (26) provides in a way a 
substitute to the probabilistic assumption of "spatial homogeneity of the flow" (see Theorem 6 and §6.2) 
and will also prove to be well suited to unbounded domains (see §3, Theorem 9). 

Examples 

1. If is bounded, the Cauchy-Schwarz inequality provides for any / 6 T 2 ([To,Ti] x Q : 

Vol(/; [T , Tl]) < Vol(fi). 

2. If f(t,x) = tn>+r)(t)( x ) is the characteristic function of a subset W of fi translated by a vector rj(t) 
such that n' + rj(t) C ft, then one has Vol(/; [T ,Ti]) = Vol(fi') and 

v„ 1(/ + £(1 - / , ;l r ,T ll) .((l±^)v„ W ^ ,= M. 

3. A simple computation on M 3 gives : 

Vol(e^ ; [To.Tx]) = Sy^vr 3 / 2 j^J^L . (27 ) 

In particular the volume is on [0,T] because 5o £ 2 and infinite on [T, oo]. On [T, AT] it is of 
the form C\(vT) 3 / 2 which conforms to the intuition that the heat kernel is around time T, mostly 
concentrated in a sphere of radius 

4. The velocity flow associated with an inviscid vortex line ui = 8 x= q® 5 v= q®1-\ <z< \ behaves as , w 

V 1 +v 

along (0,0) x [—1,1] hence belongs to L 1 ([— 1,1] 3 ) but is not square-integrable on [— 1,1] 3 ; one has 
therefore Vol^^ ^p) = 0. 
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5. The infinite viscous Oseen vortex line of direction e3 is the solution of (1) given by 

u(t,x) = (u h (t,x),0) G M 3 

with 

u h {t,z) = ±v el 2 and «(0 = f- |^(l-^ |2/4 ), reE. 

The vorticity is the 2D heat kernel u?(i) = (4irvt)~ 1 e~( x i +x ^/ 4 ' ut e%. Its characteristic scale is Ay/vi. 
As it is constant along the z-axis, it does not belong to any L P (M 3 ). One can however easily compute 
the volume function in restriction to the cylinder VL = {(x,y,z) ; x + y 2 < 1, \z\ < 1}. The result is 
shown in Figure 3. 

At the scale of the vortex still appears concentrated around the z-axis at t = 5 x 10~ 3 . The peak of 
the volume function around t = 7 x 10 -2 occurs when the characteristic scale of the vorticity matches 
that of O. 

Conversly, this simulation illustrates that for a given t > 0, the length A = Vol(w : [t/2, t]) 1 / 3 (computed 
on a large enough domain) determines the characteristic scale Ayfui of the vorticity. Numerically, one 
has 

Vol(n:[t/2,t]) / 4y^t \ 1/4 

Vol(fi) ' VVol(n)V3 J 

in the decades before the volume function reaches its peak (i.e. in physical terms, when the character- 
istic scale of observation Vol(O) 1 ' 3 exceeds the characteristic scale of the vorticity) . 

Voi(u)/Voi(n> 




Figure 3: Plot of ||«(£)||x,i / 11^(^)11x2 and of Vol(u : [t/2,t]) for V = v = 1. The computation is done in 
restriction to the cylinder = {(x,y,z) ; x 2 + y 2 < 1, < z < 1} and the result is displayed in Log-Log scale 
as a percentage of Vol(O). Inlaid pictures represent the vector field Uh(t,x) for t G {5 x 10~ 3 ,7 x 10 _2 ,10 2 }. 

1.6 Turbulence in the spectral sense of Obukhoff-Kolmogorov. 

This section provides the mathematical background of Obukoff-Kolmogorov's spectral theory of tur- 
bulence, known as "K41 theory" in reference to the publication date of Kolmogorov's [40], [41], [42] and 
Obukhoff's [57] founding papers (see [31, p. 98] for a precise chronology). 
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1.6.1 K41-functions 

Let us start with an abstract mathematical definition. 

Definition A function u G £(fi) = L°°(M + ; L 2 (Sl))nL 2 (M+; H 1 ^)) is said to be a K41-function on [T , T-y] 
if there exists C g]1, 2[ such that 



[ K 2 E\K) d^{K) < C [ K 2 E\K)dn{K) 



(28) 



where E(«; [T ,Ti]) = G ^(A 1 / 2 ) ; K > Vol(u; [T , Ti])" 1 / 3 } and Vol(u; [T ,Ti]) is defined by (26). 
On M 3 , if n is a K41-function then for any K_ > such that 

(K_) 3 x Vol(«;[T ,ri]) < 1 (29a) 

one has for all K + > K_ : 

*K 4 



POO r i v _|_ 

/ K 2 E*(K)dK < (l + C) K 2 E*(K) dK. 

JO JK_ 



(29b) 



Any such interval [ET_,i^ + ] is called an inertial range of u. The corresponding spectral Reynolds 
number is : 



5ft 



(29c) 



One also defines a spectral precision parameter (the reason of the fraction 5/3 appears in the next section) : 



7 = sup 

Ke[K-,K+] 



(29d) 



Obviously 5ft and 7 are dimensionless and there are infinitely many admissible quadruplets (K±, 5ft, 7). 
On T 3 , the smallest non- vanishing frequency possible is 2nL~ l so (29a) is replaced by 

Vol(u;[To,Ti]) 



X_ = 2vrL^ 



if 



Vol(T 3 ) 

2nL- 1 <K < Vol(u; [T , T1])- 1 / 3 otherwise. 



> (2vr)- 3 , 



(30a) 



Note that one always has Vol(u; [To, Ti]) < Vol(T 3 ) by the Cauchy-Schwarz inequality. The constant is 
(27r)~ 3 ~ 4 x 1CT 3 . In particular, one has : 



2vr < LK_ < max { 2vr; 



Vol(T 3 ) 



Vol(u;[r ,ri]) 



1/3' 



Admissible frequencies K + are defined by (29b) with an obvious change of notations. The discrete substitute 
for the definition of the spectral precision is : 



7 = max 

K-<K j <K j+1 <K + 



i0 § E^Kj) 5 



(30b) 



where E = a(A 1 / 2 ) = {Kj ; j G N} with /C,- < K j+1 . 
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1.6.2 K41-turbulent flows 



Let us now turn back to fluid dynamics. 

Definition Turbulence in the Kolomogorov-Obukov sense is the question to find and describe solutions 
of (1) that are K41-functions on some time interval [To, T±] and that possess at least one inertial range K± 
in the asymptotic regime : 

K>1 and 7 logSR<l. (31) 
Such a solution is called a K41-turbulent Bow. 

In section §2, the asymptotic (31) will be used to recover the K~ 5 / 3 law found in physics textbooks. 
Remarks 

1. It is a mathematically open problem to construct exact solutions of (1) that possess this definite 
behavior. However (31) has been observed in numerous experiments as the generic state of highly 
fluctuating flows (see e.g. [8], [62], [7], [53], [67] and the references therein). Physics textbooks usually 
don't mention the spectral precision 7 because they rely on Log-Log plots of the energy spectrum on 
which the property 7logK <C 1 is equivalent to having a substantial amount of data concentrating 
along a straight line of slope —5/3 on logK decades of frequencies ; 7 is the relative error on the slope 
of the line. 

2. In naive terms (29b) means that, by definition, the K41-theory of turbulence is a spectral property 
of vortex structures because at least half of the average enstrophy (IMI^) comes from the frequency 
range [K-,K + ]. 

3. The use of Vol(u; [To, TJ) at the foundation of K41-turbulence is new. There is experimental evidence 
that K41-turbulence is generated on some thin-structured subset of 0, deeply connected with the 
vorticity and whose characteristic size £q = K- is the largest scale involved in the inertial range. 
Example 5 p. 11 has already shown a strong but subtle connection between the volume function and 
the characteristic scale of vortex structures. It will be shown (see Theorem 6) that 

^Vol^fTcT!]) 1 / 3 

if u is a K41-turbulent solution of (1). 

Let us conclude this section with a short dictionary between mathematics and physics. In physics 
textbooks, one can find statements like 

<C Turbulent solutions u of (1) satisfy F(u) < G{u) > (32) 

where F and G are two functionals on C{Q) possibly depending on To and T\. From a mathematical point 
of view, it should be read as follows : there exists a function C : R 2 — >]0, 00 [ with 

lim C(»,7) = C €]0,oo[ 

5R-5>oo 
»t-»-1 

and constants 3?o,7o > such that any solution of (1) that is a K41-function on [To,Ti] with a Reynolds 
number !R > 3?o and a spectral precision 7 < satisfies 

F{u) < C(K,7)G(«). 

In this case, any solution that admits parameter (3?, 7) in the asymptotic range (31) will indeed satisfy 

F(u) < C' G(u) with C' /C ~ 1. 

One should be aware that that physics textbooks usually contain additional "meta-assumptions" such as 
the isotropy or the homogeneity of the flow, which should then be translated adequately and added to the 
assumptions of the mathematical statement. 
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2 General properties of K41-functions 

2.1 Scale invariance and stability in Leray space 

The definition of K41-functions is invariant under (20) : If u is a K41-function on [To,Ti] with parame- 
ters (if±,§R,7) then u\ is a K41-function on [A" 2 To, A" 2 Ti] with parameters (A-?T±, 3ft, 7). 

The set of K41-functions is open in the Leray space £(T 3 ). More precisely, the following statement holds. 

Proposition 4 Assume that u G £(T 3 ) is a K41-function on [Tq,Ti] with parameters (K±,^R). For any 
e > there exists C e > such that any v G £(T 3 ) with 

ll n ~ v IIl 2 ([T ,Ti]xT3) + II ~ t, )lli 2 ([T ,T 1 ]xT3) ^ ^ 

is aiso a K41-function that satisfies (29a) and (29b) with the same parameters (K±,$t) but with numerical 
constants 1 + e and 1 + C + e. 

Proof. Each term of (29a) and (29b) is continuous on ^([I^Ti] x T 3 ). ■ 

Note however that the spectral precision 7 is not preserved, which means that if a flow u satisfies (31), 
any neighborhood of u in C(fl) will also contain functions that do not satisfy this asymptotic. This illustrates 
a wise comment by U. Frish [31, pp. 199-202]: "Questions (on turbulence) are likely to benefit from a close 
collaboration between mathematicians and physicists but it will require more than better functional analysis 
(...) ; some geometry is needed." 

2.2 Kolmogorov's constant a and the K~ 5 / 3 law 

Property (29b) involves the dissipation rate e (which is the left-hand side) and frequencies. There is only 
one way to define a quantity that has the same units as the energy spectrum [length] 3 • [time]" 2 (i.e that 
scales the same way under (20)) and which is a power function of a dissipation rate e = [length] 2 • [time] -3 
and of a frequency K = [length]" 1 : 



^2/3^-5/3 = [ length ]3 . [time j 



-2 



This fact makes the | fraction in (29d) a more obvious choice. Note that even though E*(K,t) and e{t) 
depend on p, the dimensionless fraction 

E*(K,t) 

e (t)2/3#-5/3 

is not "missing" a power p 1 / 3 because it is invariant under the transformation (20), which does not change 
the total mass. 

The following property of K41-functions is often mistaken for a definition of turbulence. Its real meaning 
according to the dictionary (32) is that the most valuable theorems concerning K41-functions will hold in 
the asymptotic regime (31). 

Proposition 5 If u is a K41-function on [Tq,Ti], the function a(K) = _ 2/3 ^_g /3 satisfies 

VK, K' G [K- , K+] , K" 3 ^/ 4 < ai ^X < » 37/4 (33) 

a(K') 

Proof. The definition of 7 (namely (29d) for M 3 and (30b) for T 3 ) implies : 

E*(K\) (K^ 
K ) ~ E*(K ) ~ \K 
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for any Kq < K\ in [K-,K + ]. Applying this inequality either to (K,K') = (Kq,Ki) or to (Ki,Kq) one 
gets : 

\K-J ~ \K>) ~ a(K>) ~ \K'J ~ \K-J 
hence the result. ■ 



Definition In the asymptotic regime (31), the function a is essentially constant on the inertial range and 
is called the Kolmogorov constant. To fix further computations, one will choose from now on : 

a = a(K + ) (34) 

±1 



and (33) reads ( ^3^5/3 ) <^ !i on[K-,K + \- 



Remark Proposition 5 would also hold if in the definition of ot(K) one would replace e by another quantity 
having the dimension of a dissipation rate, which in turn would change the value of the Kolmogorov con- 
stant (34). There is physical evidence that this definition leads to a universal numerical value for a (see §6.5) 
but mathematicians should question it. The author is grateful to W. Craig for this valuable remark. 

2.3 Bounds of the inertial range - Expression of K± 

Using dimensional analysis, there is only one way to define a frequency as a function of e and u : 

1/4 



K d = a" 3 / 4 



U • < 35 > 



Since a is dimensionless, the power of a is arbitrary here, but it is the one that provides the simplest state- 
ment in the following Theorem 6. It also corresponds to the "phantom" homogeneity a ~ p 1 / 3 mentioned 
above. In physics, the length n = K^ 1 is often referred to as the "Kolmogorov dissipation scale". 
Likewise, a frequency defined as a function of e and E is : 

Kc = a^^j- 2 . (36a) 

In physics, the length £q = K^ 1 is refered to as "the size of large eddies". For T 3 , the energy spectrum 
misses the total impulsion (see (12) and §3) so one substitutes the following definition for (36a) : 



K c = a 3/2 e [ E 



puo(x)dx 

T 3 



2^ 



-3/2 



(36b) 



Physics textbooks usually state that K + = and K_ = K c by computing E and e for an idealized 
compactly supported energy spectrum on M 3 : 

E*(K) = ae 2 / 3 K' 5 / 3 l [K _ tK+] (K) 

Converting this idea into a rigorous proof must be done carefully and actually requires some additional 
assumptions to hold for K_ . Moreover, the computation on T 3 (which always seems to be dodged in physics 
textbooks) happens to be extremely instructive (see also Remarks 1 and 2 p. 6). 

Theorem 6 The following inequalities hold. 
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1. Case of SI = M. 3 . Ifu is a K41-function on ft = R 3 with parameters (K±,$l, 7), then : 

( 3?~ 37 / 4 V /4 K+ ( 4K 37 / 4 V /4 

< IT < on SZri - (37a) 



6(1 - K- 1 ) y - K d - \ 3(1 - K- 1 ) 
3/2 ^ 



Moreover, if u is a solution of (1) with initial data uq € L 1 n L 2 (M 3 ) then 



12 SR3 7 /2 - V -y v .L u, u/_ I , V J Vol ^3 

Corollary 7 In the turbulent asymptotic 3? — >■ 00 and 5ft 7 — > 1, one has 



for two numerical constants C, C . 



(37b) 



^^(^x^l-^V 72 (37c) 



and 

!:i ' T " ~ hK " 2 < (^-) 3 x Vol(u;[To,Ti]) < 1. (37d) 



36(1 - 5? -1 / 2 ) 

2. Case of SI = T 3 . Ifu is a K41-function on £7 = T 3 with parameters (K±,$l, 7), then : 

<^< 1, , 1 — > W 

V 6 J ~ K * ~ [ if - 1^ 1 J 

( K 37/4 (if " 35rl/2 )) 7 - ^ provided K- > S^irL- 1 ). (38b) 

Moreover, if u is a solution of (1) with / puo{x)dx = and 

„, , Card{zGZ 3 ; |z| 2 = n_} f LK -Y /on , 

C(n_) = 1 < 1 where n_ = , (39) 

v ; 8vr 3 ^/nT V 2vr / 

then : 

K_ / 12K 37 / 4 V /2 , . 

and 

1 " C{U - ] < (K-f x Vol(«; [To, rj) < max (l; (2vr) 3 . (40b) 



Vol(u; [T ,Tij) x /3 



Assumption C(re_) < 1 is still an open problem in number theory ; the numerical test presented in 
Figure 4 ensures it is satisfied for most if not all practical purposes. 
Numerically, the theorem reads for Q = M 3 : 

0.260.Kd <K + < 1.25K d , l.83K c <K< 6A6K C and 0.711 < {K_f x Vol(«; [T , Ti]) < 1. 

For VL = T 3 , the asymptotic (31) implies K + — > 00 because if_ > 2itL~ 1 , thus (38a) provides — > 00. 
One has C(n_) < 0.1 for n_ < 10 5 and 

0.260^ <K + < 0.625JQ, 0.907if c < K_ < 48.7i-C c 

and 

8.33 x 10~ 3 < (K_) 3 x Vol(u; [7b, Ti]) < 248.1. 
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Remark One could object that the numerical constants of the previous statement are not fundamen- 
tal. Indeed, for most of what will follow, one could just write C(JR, 7) with C(Jft,j) — > C > in the 
asymptotic (31). However, for Theorem 16, which investigates the subtle relation between smoothness and 
intermittency, the numerical values of the constants will mark the difference between an empty statement 
and a meaningful result, which means that one will have to show some discipline in each intermediary result. 

Proof (estimate of K + ). One can compute e using (17). On the inertial range (33) leads to : 

S > „ £ K 2 E*{K)dK >_ ^(tff - K<f>) = - ^tff 

hence the right-hand side of (37a). Conversely, according to (29b) one can compute e using only the inertial 
range where (33) again provides : 

r-oo r-K + 

e<Au K 2 E*(K) dK < 8crf 7/4 e 2/3 / K 2 ^ 3 dK = 6crf 7/4 e 2/3 (l - ST 1 )^ 3 
Jo JK- 

hence the left-hand side. On T 3 , one gets instead : 

KeT,f][K^,K+] ^ ' 

hence 

K -3 7 /4 fl < (*A 2 y K -2/3 < til tfhl\ 

4ai> V L J ^ lav 

The next step is to show that the center term is equivalent to KJ (the apparently different game of powers 
reflects the spectral asymptotic (14) of the Stokes operator ; note that physical dimensions are the same). 



rard|/fcy';|/| ; = n Vn") 




Figure 4: Computation of C(n). 

Numerical test that r 3 (n)/(87r 3 v / n) < 0.1 for n < 10 5 
where rs(n) = Cardjz 6 I? ; \z\ 2 = n}. A similar test 
for 10 10 <n< 10 10 + 100 checks the same numerical 
bound C(n) < 0.1 with a maximum of 9.37 x 10~ 2 for 
n = 10 10 + 1. It is obtained using SquaresR[3, n] 
with Mathematics©. The problem of computing the 
number of representations of an integer as the sum of 
three squares has been addressed historically e.g. in [27] 



and [3] but the asymptotic behavior of r^(n) cannot be 
read directly on Hardy's explicit formula. The series 
En<7V r 3(n) = fiV 3 / 2 + 0(7V e + 29 / 44 ) is the num- 
ber of lattice points inside the sphere of radius y/~N 
(see [16]). Similarly, a recent paper [21] shows that 

£„<W™) 2 = Ms) n2 + °( Nl4/9 ) but a s ain tne 

remainder is too large to prove the asymptotic behavior 
^{n) = 0(y/n). Note that this conjectured asymptotic 
is extremely sharp because the sequence r^n) vanishes 
on the subsequence r^(A p (8q + 7)) = thus cannot 
have a power law equivalent. The best known estimate 
[16] is r^(n)/y^n = 0(n e ) for any e > 0. Thus it might 
be that r^(n)/^/n is not bounded but grows extremely 
slowly. However, from our numerical data, one can infer 
that even if it diverges as logn, then r^(n) will finally 
exceed &ir 3 ^/n only for n ~ 10 50 . Such a scale is unreal- 
istic because it would require the building of a periodic 
torus of at least L = 6 light years in order to ensure 
that such a large wave number would investigate scales 
that exceed the atomic one (27rL _1 y / n)~ 1 = 10 _9 m. 
Therefore, in the range of validity of Navier-Stokes 
equations, assumption (39) is numerically satisfied. 
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-1/3 



One has : 

(?) 2 E — (t) 4 ' 3 E - 

v 7 A"esn[.K"_,.K" + ] v 7 rien 3 n[n_,n + ] 

where ^3 = N\{4 P (8(/ + 7) ; p, q £ N} denotes the set of integers that are the sum of three squares. Its 
complimentary is included in the subset of integers n = 0, 4 or 7 mod 8, therefore : 

N N 

n=l nen 3 n[l,JV] ra=l je{0,4,7}8fc+j<JV 

which after comparison to an integral boils down to 

3Nm >- E ^,|(a + ^-.)-^(!M)-^ + ow 



2 

nen 3 n[i,AT] 



Af-5>oo 



V 2 / 

One uses those estimates to compute the sum for K G E n [.ftE , : 

Note that in the asymptotic (31), one has K + — > 00 because K- > 27rL~ 1 . One finally gets : 
3T 3 ^ 4 ^ < -< 3 and < 3 ( £j - I®' 1 ) + O(l) < ^ K 37//4 . 



Proof (estimate of K-). To estimate K-, one compares the total energy with the energy contained in 
the inertial range. If u is a K41-function on M 3 , one has : 



poo poo /*-"■+ 

/ E*(K,t)dK < R- 2 K 2 E*(K,t)dK < R- 2 K 2 E*(K,t)dK < E*(K, t)dK 

JK+ Jk+ JK- JK- 

and (3b) provides 

pK+ f'K+ pK— 

/ E*(K)dK<E<2 E*(K)dK+ E*(K)dK. (42a) 

JK- JK- JO 



On T 3 the corresponding inequalities are 



fl"esn[.K"_,.K" + ] 

2 



* - (X, E (il) fi * (A ' >+2 E (£) *w 



(42b) 



Using (33) on [K-,K + ], the lower bound of E provides (37b) : 

,2/3 3^2/3 



The computation is similar on T 3 ; one gets 

2/3 



_zj > sft-37/4 x n 1 J 3 c(n-,n+) with C(m,n 2 )= 



neD3n[ni,n 2 ] 
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This time, one uses C(n_, n + ) = C(n_, oo) — C(l + n+, oo) with 



3 (AT _ l)-l/3 > C(jV) > ^AT-1/3 + 0(^-4/3) > 



the last inequality being valid if N > 10. The proof of the converse inequality and the last statement about 
Vol (it; [To,Ti]) will be postponed until Proposition 11. One will use the additional assumptions and (29a)- 
(30a) to estimate the spectrum for low- frequencies in (42a)-(42b). ■ 

3 Low-frequency spectrum 

In physics, two different low-frequency spectra have been described (see e.g. [25, p. 91 and p. 358]) : 
the Saffman spectrum behaves like K 2 and the Batchelor spectrum behaves like K 4 as K — > 0. Numerical 
simulations on T 3 have confirmed that both can exist and that the choice between K 2 or K 4 depends on the 
initial data : the initial behavior of the spectrum is preserved, though in case of K 4 the spectrum does not 
seem totally stable in time (see [51, fig. 5a and fig. 7] and the references therein where this subtle instability 
is referred to as backscatter). For real- world experiments, the question is not properly answered because 
the first Fourier coefficients reflect the large scales in the experimental protocol more than turbulence itself. 

In this section are proved the following points : 

• Theorem 9 : Saffman's estimate holds in general. 

• Theorem 10 : On R 3 , the spatial localization of the initial data determines spectra in K 2 ^ l+ ^ with 
possibly any (3 < 1 depending on the exact decay at infinity. 

• Batchelor's K 4 spectrum (/3 = 1) seems to occur only for unstable highly localized solutions. Fur- 
thermore, it was not possible to obtain any nontrivial estimate with f3 > 1 even for highly localized 
flows. 

Note that on T 3 , one can always artificially improve the power of K £ X* of upper bounds by multiplying 
by KL/(2ir) > 1, thus it is worth pointing out that (45) below does not contain the length L on the 
right-hand side. 

3.1 Estimate from u e L 1 (fi) 

From a physical point of view, it is reasonable to assume that 



which means that the fluid is globally at rest in the given coordinate system. On Q, = T 3 , Cauchy-Schwarz 
ensures that u(t) G L 1 (il) for any t > ; moreover, the total impulsion is constant in time so (43) holds for 
any t > if and only if it holds for the initial data. On M 3 , the following statement justifies why (43) holds 
anyway. 

Theorem 8 If u is a Leray solution of (1) with initial data uq G L 1 n L 2 (R 3 ), then u(t) 6 L l (Q) for any 
t > and (43) holds for (at least) almost every t > 0. 

Proof. On M 3 , the following inequality holds for Leray solutions of (1) : 




(43) 





(44) 
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which ensures that u(t) G L 1 (M 3 ). Let us conclude first before proceeding to the proof of (44). Leray 
solutions are known to be smooth for almost every t > (see [50]). At any such time, the following 
integration by part holds for any smooth bounded subset SJcl 3 : 

/ Ui = — Xi(divu)dx+ / XiU-dn. 
Jn Jn Jan 

For such times, the fact that u(t) G L 1 n C°(K 3 ) also implies that : 

I I Xiu(x) ■ dn(x) 

>JR J\x\=r 



lim 



*=0. 

r 



One can choose a sequence of radii — > oo such that lim / X{ u(x) ■ dn(x) = 0. Combining both 

*->-oo 7|x|=r fc 

identities and divu = gives (43) for any time t > such that u(t) G H 1 n C°. 

Inequality (44) is part of the folklore of fluid mechanics and, contrary to most estimates, it does not 
involve an exponential growth in time. Let us recall briefly its derivation. Rewriting (1) as a heat-equation 
with a non-linear source term (see [48, chap. 11]), one gets : 

u(t,x) = e utA u (x) + f [ (u®u){t',x')K(t -t',x - x')dx dt' 
Jo JR 3 

with a convolution kernel that satisfies \K(t,x)\ < C(\x\ (see e.g. [66]). Therefore, for any r G [0,t], 

one has : 

ft f°° r 2 drdt 

\\u(t)\\ l1 < \\uo\\ L i + c||u®u|| L oo ([ o it ] ;£ i) y o y o 

|2 



/ (r + Vz4) 4 



But 1 1 it ® w|| £ i < 1 1 wo 1 1 1,2 and the last integral computes down to hence (44). ■ 

One can now deduce bounds on the lower end of the energy spectrum. 

Theorem 9 If u is a Leray solution of (1) with u G L 2 n L 1 ^) on Q = R 3 or T 3 , then for any T < T\, 
the energy spectrum on [To, T\] satisfies : 

E*(K) < Vol(«; [T , Ti])K 2 E x 4/vr 2 on R 3 , 

(45) 

E*(K) < Vol(u; [T , Ti^E x Card{z G Z 3 ; \z\ 2 = n} on T 3 , 

where the volume Vol(u; [T ,Ti]) is defined by (26). For T 3 , one has L = Vo^T 3 ) 1 / 3 and K = 
with an integer n G CI3 that is the sum of three squares. 

In the case of fi = T 3 , one does not require (43) to hold. The numerical factor is illustrated by Figure 4. 
Proof. Using ||n|| ioo < ||w|| L i, Proposition 1 reads on R 3 : 

1 r°° 

E*(K,t) <p\\u(t)\\ 2 Ll K 2 x ^r)r 2 dr 
The integral is bounded by 8 because suppV' C [2 _1 , 2] and f ijj(r)dr/r = 1. On T 3 , one has instead : 

E^K,t)<p\\u(t)\\ 2 Ll K 2 x ^-^L. Card{fc G f 3 ; \k\ = K}. 
Conclusion (45) then follows immediately from the definition (26) of the volume function. ■ 
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3.2 Estimates from localization norms of u on IR 3 

One can improve the low-frequency estimation of the spectrum using localization properties of the flow. 

Theorem 10 Given (3 g]0, 1[, there exists a constant Cp > such that any Leray solution of (1) on = M 3 
with initial data in L 2 n L 1 satisfies for almost every t > : 



L 2 (R3) 



(46) 



(47) 



E*(t, K) < pC p K 2{ * 1+ ^ (1 + \x\f 2+ P U {t) 

If the initial data uq satishes e.g. 

(1 + \x\f +l3 'u G L°°(M 3 ) with /3'>/3 
then the right-hand side of (46) is hnite as long as u is a smooth solution of (1) on [0, t]. 

Remarks. 

1. For (3 = 1, the estimate (46) holds with an extra multiplicative factor (1 — log.fr/-Ko) for K < Kq but 
the right-hand side is infinite unless a non-generic (necessary but not sufficient) condition holds : 



/ 



2 

Ui(t,x)uj(t,x)dx = - |Ki)||x,2 Si j. 



(48) 



This condition is not invariant by the flow of (1) and [13] contains examples of flows that will check 
and violate (48) at prescribed times. Those smooth flows will satisfy (46) for any f3 < 1 and t > but 
they will satisfy it for /3 = 1 (with the logarithmic correction) if and only if t G {to, t±, . . .}. Examples 
of flows that satisfy (46) for all time with (3 = 1 (with again the logarithmic correction) can be found 
in [11]. 

2. Contrary to (44), the best known short-time bound for weighted norms is exponential in t. Other 
localization norms of uq could be used instead of (47) . For example the same theorem holds if : 



(1 + lxiyH 1 p) +p 'u (x) G L p (R 3 ) 

with p > 3 and (3' > /3 (see [66]). 

3. On £2 = T 3 , weighted norms are meaningless ; however, provided (43) is satisfied, one could get a 
similar result with the weighted norms replaced by : 

\u(t,k) - u{t,k')\ 



I«(*)llcw = SU P 
k^k' 



\k-hf\P 



To the best of my knowledge, propagation of this semi-norm by the flow on T 3 for (3 g]0, 1[ has not 
yet been studied. Neither has the propagation of u(t) G J 7 ~ 1 (C /3 (IR 3 )) in the continuous case. One can 
however expect a generic failure of the propagation of sup^g |A;| _1 |n(t, k)\ corresponding to /3 = 1. 

Proof. One has only to prove (46). The rest of the statement follows from [66] for (3 g]0, 1[. One has the 
following chain of continuous inclusions if f3 £ N : 



\u\\ C f> ^ IHIjyl+iS < 



{l + \x\f /2+f3 u{t) 



L 2 



If p G]0,1[, one gets \u{t,£)\ < C^l 13 ||(1 + |x|) 3 / 2+/3 u(t)|| L2 when n(t,0) = j R3 u(t,x)dx = 0, which (accord- 
ing to Theorem 8) holds for almost every t > 0. Therefore : 



E*(K,t) = (2tt)- 3 pK 2 / / 

Jo Js 2 



i>(r)\u(t,Kr$)Yr 2 drdf) 



< 



CpK 2 ^ (jH \p(r)r 2 ^dr^j || (1 + |x|) 3 / 2+ ^(t) 



2 

L 2 
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which gives (46) with a constant independent of ip because supp^ C [2 1 ,2] and f ip(r)dr/r = 1. 

Let us briefly justify the notes that follow the Theorem. The precise relation between (48) and the 
finiteness of the right-hand side of (46) is extensively studied in [12]. When (3 = 1, the Sobolev space 
H 5 / 2 (R 3 ) is not included in Lip(M 3 ) but in Calderon's space Cl(R 3 ) for which - f(rj)\ < C|£ - fj|(l - 
l°g(|£ ~~ ^D) when |£ — 17 1 < 1 ; see e.g. [17, p. 31]. When (3 > 1, e.g. for /3 e]l,2[, the following estimate 
holds : 

\u(t,0\ = \u(t,O-u(t,0)\ < |£| \Vu(t,aO\da^J < \£\ ^f' 1 \\u(t)\\ c , j\^da + \\\x\u(t)\\ L ^j 

so if u is a highly localized flow like those in [11], then (46) holds with (3 > 1 ; however, K 2 ^ +1 ^ must be 
replaced by K A because |u(t,£)| / o(|£|), even for highly localized flows. ■ 



3.3 Precisions on K_ 

Low- frequency bounds on the energy spectrum allow us to conclude the proof of Theorem 6. 

Proposition 11 Let u be a Leray solution of (1) and a K41-function on [Tq,Ti] x fi. 
1. IfQ = R 3 and u G L 1 n L 2 (M 3 ), one has 

t £ (^ xS3T/1(1 - s " 1/2) ) 3/2 (49) 

and 

:5:r '"" : <(# ) :, x Yol(, /:"/;,. V,V 1. (50) 



365R 3 T/2(l-SR-l/2) 

2. Ifft = T 3 , let us assume that f pu (x)dx = and 

Jt$ 



Card{z G Z 3 ; \z\ 2 = n_} f LK ^ 2 



C(n_) = „ < 1 where n_ = 



•S-\ /; V 2,7 

Then similar inequalities to (49) and (50) hold, namely 

K-^y^wi) (51) 

and 

1-C(n_) . ^3^ 1 , - ^n^^/i./o^sVol^tTcT!])- 



12K 3 7/2 



< {K_f x VaL(«; [To,^]) < maxjl; (2vr) 3 ^/ff" } • (52) 



Assumption C(n_) < 1 is still an open problem in number theory but the systematic numerical test presented 
in Figure 4 ensures that it is satisfied for at least n_ < 10 5 , which should be sufficient for most practical 
purposes. For example, in a periodic domain of size L = lm, it means that the assumption is satisfied at 
least as soon as the size of large eddies KZ 1 exceeds 0.5mm. Note that there are no restrictions on K + , 
which means that turbulent structures can develop details at much finer scales. 
Proof. On M 3 , the starting point is (42a). Using (33) on [K-,K + ], one has : 

J B<3a£ 2/3 (^_)- 2/3 K 37/4 (l-K- 1/2 )+ / E*{K)dK. 

Jo 

Then (45) gives E*{K) < ^VEK 2 with V = Vol(u; [T ,Ti]), hence : 

/ E*(K)dK < —^VE(K-) . 

J0 
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Then (29a) reads V{K-) 3 < 1, so £ < Zae 2 / 3 {K-)- 2 / 3 ^/ A {l-^ l / 2 ) + ^E and (49) follows immediately. 
To get (50), one checks the compatibility between (33) and (45) at K = K_, which provides : 

$l- 3 ^ 4 ae 2 / 3 {K_)- 5 / 3 < E*{K_) < 4j Vol(u; [T , T^KJfE. 

IT 

On recognizes ae 2 ! 3 = EK^ 3 , hence using (49) : 

x Vol(„; [To.TJ, > ^ > ^ x - ^)) . 

The upper bound was already given by (29a). 

For T 3 , the starting point is (42b). As before one uses (45) on [27rL _1 , K-\ : 

3 



(¥) E (^) ^ ■ (t) \£^y*MTo,Ti])E 

< C(n_) • (x^ 7 ) Vol(u; [ T °' Tl ])^ 



with the exact same constant 
= C(n_) • (K-) 3 Vol(«; [To,Ti])S. 

If V ° > (2tt) _3 then (30a) gives -fT_ = 27rL _1 and the sum on the left-hand side is empty so there is 

nothing to estimate. If not, then (30a) gives (K_) 3 Vol(u; [To,Ti]) < 1 and the sum is bounded by C(n-)E. 
Since it is assumed that C(n_) < 1, one can bootstrap this term in the left-hand side. 
On [K_,K + ] one uses (31), which gives as before : 



n[K-,K + ] v 7 I nen 3 n[n_,n + ] 



< 6K 37/4 ae 2/3 KI 2/3 if n_ > 2. 

For the last inequality, one has just estimated the sum on n [n_, n + ] by the sum on all integers greater 
than n- and then compared it to an integral. Then (51) follows from : 



(1-C(n-))E< puo(x)cbj + I2K 37 / 4 | E - (^j pu (x)dx 



2 ) K-- 2 / 3 

( , ptiiA.r \<l;v ! 



Compatibility between (33) and (45) at K = K_ then gives : 
using again the assumption C(n_) < 1. 



4 High-frequency spectrum. 

Physics textbooks often state that turbulent spectra have rapid decay at high-frequencies, i.e. that 

ViV G N, sup K N E*(K) < 00. (53) 

K>0 

Property (53) implies that u{x) is smooth. Conversely, the following result relies on the best known smooth- 
ing effect for (1) and implies that (53) is automatically satisfied for smooth solutions. 
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Theorem 12 There exist (dimensionless) numerical constants C, Co > and a family c(K) > with 

' J2 C ( K ) = C if = T 3 

KeT,* 

\/K G R, ^ c(2 n if) < C if = R 3 

that iiave tiie foiiowing properties. For any smooth solution u of (1) on [Tq,Ti] x Q, with u(Tq) £ 
let us define 

and S(t) = - min{ y/v(t-T ); C y 7 ^} . (54) 



su P[T ,Ti] l|Vu||£ 2(Q) ' ' 2 

Then for any t 6 [To, Ti] the energy spectrum satisfies : 
• on SI = T 3 : 

G £*, £ f (K,t) < C (K) • ^Vol(T 3 )^ e-^Eo, 



(55a) 



• on n = R s : 

r2K 



i r 

\/K et;, — / E*(t, k)dk < c(K) • K~ x e- &i ^ K E , (55b) 
with Eq = E(Tq) the initial kinetic energy. 

Exponential decay in (55) with a uniform constant 5 on [To + r, T\] means that the analyticity radius of 
smooth solutions remains uniformly bounded from below. This question has raised some concerns in the 
asymptotic v — > (see [31, p. 92-93] and the references therein for a brief survey). For the question of 
analyticity for a given v > 0, the end of this section contains detailed bibliographical notes. 

Even though the main interest of the statement is the high-frequency asymptotic, it should be compared 
to the low- frequency inequality (45) on T 3 . 

Corollary 13 If u is a smooth solution of (1) on [T , T{\ x T 3 then for all K £ S* : 

&{K) < C ( ^ Vol(T 3 ) J E <^- Vol(T 3 )K 2 T . (56) 
\ L J 2tt 

Proof of Theorem 12. The idea is to first prove a short-time analyticity estimate (61) valid for any 
Leray solution of (1) evolving from smooth initial data. Then one uses the decay of kinetic energy E(t) and 
the qualitative assumption that u remains smooth on [0, A] to propagate this quantitative estimate along 
the time line. 

Let us consider the following ODE with unknown function fl(t) and a dimensionless constant A that will 
be adjusted later on : 



0(t)= [ \C\e Mt ' m u(t',C) 2 dt' with A(t) = V^i- 



Ad(t) 



(57) 



To be perfectly rigorous, one should consider a family of smooth approximations u n of u, e.g. Friedrichs 
approximation by low-pass filters. Then (1) and (57) are of Cauchy-Lipschitz type and can be solved for all 
t > 0. The auxiliary ODE (57) will be used to prove an analytic estimate (61) of u n that will pass to the 
limit u n — > u on some uniform time interval. To keep the formulas reasonably sober, the index n is dropped. 

The image of (1) by the Leray-Hopf projector P is 

dtu — uAu + P div(n (g) u) = 0, u(0, x) = uq(x). 



24 



DuhamePs formula reads 



and 



Mt)\C\, 



u(t,0 <W L (t,0+ I W NL (t,t',Odt' with 
J o 



W L (t,0 = e 



w NL (t,t',o = |^| e [ A w- A (*')]i«i-^')i€i 2 ^ e Mt>m\ J ^ lu^^-r?)!!^',!?)!^ . 

To simplify notations of the time integral, one defines also : 

Wnl(*,0 = I V nl (m' ,6^' '• 

JO 

The same formula holds on T 3 with the obvious modifications in the notations of discrete spectra. Let us 
first estimate the linear term. One has for any t > : 

Ami - »m 2 < v^tiei - ^iei 2 < \ - \vt\z\ 2 

therefore \\W L (t, -)\\ 2 L 2 < C \\u \\ 2 L 2 and 

= c/ (l-e-^l a )|«o(0| 2 de 

<C inf (^ 111^1^(011^ ■ 
se[o,i] 

For the phase of the non-linear term, one uses the identity (a — — ^P) < 2 : 



(58) 



[aw - A(t'm\ - u{t - oiei 2 = (jT f ) - \*{t - t> 



^\Z\(Vvt-Vvt>)[i--(Vvt + Vvt>)\t;\ 



< 



;f)-^-oiei 2 +2 . 



The sub-linearity of £ i— > |£| provides : 

e A C')l«l / \W,Z-ri)\\uy,T,)\dT,<T[V(l!,-) 2 ](Z) with V(t',0 = e A{m W,0\- 

One has llV^i)!^ = ||VV(i)||^2 = ^jf ■ A classical property of Sobolev-Besov spaces is that the product 
{fid) fd maps continuously H 1 x H 1 to H 1 / 2 ; therefore : 

o<^[^.o 2 ](o<c7^(f)-^,oiei- 1/2 with ||^,.)|| L2 <i. 

Putting everything together, the following bound holds for the non-linear term : 

\\w NL (t,0\\h<c [ /'|f|V2^-^ l J?Sw*-^-*')i€i a §(0^,0A' 

Jm.3 Jo at 



< c I (jf f(0ieie- 2 ^^ so**) (jf e"^ ICI " f (O^.O 2 *' ) « 

- A L I 2 



(59) 
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Along the same lines, one has : 



» I \\\t;\WNL(t', dt>< c " 

Jo 



< 



A 

Cv 2 '■' 



l- e -^') 



)(!> 



o 



< 



A 

Cv rt 
A 



dt 



(60) 



Estimations (58) and (60) imply the following bootstrap : 

C C 1C 

m < - M M* m\ s uo(o\\h + -jm < — m\ s ^(o\\h 

V s£[0,l] A V s6[0,l] 

provided the numerical constant A is chosen sufficiently large. To conclude the analytic estimate of u, let 
us now focus on the time interval on which A(t) > \\Jvi. Let us therefore define : 

T* = inf ji>0;Vt'e[0,i], < ^J^* | , 



The bootstrap inequality provides 



2 . 



T* > sup 



2, ,3 



> 



Civ 
l|V«o|| 



L 2 



On [0,T*], one has iy 7 ^ < A(t) < v 7 ^ hence (59) and again the bootstrap inequality gives : 

/ e^\u(t,0\ 2 d^< [ e A WI«ln(t,6 2 ^<{||^ L (t,-)IL 2 + II^NL(t,-)|| L2 } 2 <C'||no||i 2 . (61) 

This inequality proves that the radius of analyticity of u(t, •) exceeds \\fvi on [0, T*]. The constant C is 
purely numerical. This estimates also holds with an obvious change of notations on T 3 . 

Let us now assume that u is a smooth solution of (1) on [0, T] x R 3 or [0, T] x T 3 with T possibly much 
larger than T* , one can translate (61) along the time line in the following way. One defines : 

v 3 

t = j and <5o(t) = \[vt. 

sup [0 ,T] l|V«|| £ a ( n) 

Then for any time to £ [0,7"], the estimate (61) holds on [to, to + Cqt] with the same constant C : 

Vt e [0,71, V£ G [to, to + Cqt], C- 1 e §V-"(*-to)A u(v) 2 < || u (t ,.)|| L2 < IM^ . (62) 

L 2 

The best estimate at time t is therefore obtained by choosing to = t — C^r provided t > CqT. When t < CqT, 
one can only rely on the initial estimate with to = 0. 

One can now conclude the proof of Theorem 12 from a straighforward computation. As the Bessel- 
Parseval formula is slightly different let us detail it : 



'7^3 / lAOl 2 * onR 3 , 

Jr 3 



L 2 (Q) 



= < 



Vol( 



T 3 . 



(63) 



fceT s 



2 In particular, one can pass to the limit u n — > u on [0, T*] and this method guarantees smoothness. Note that T* = +00 if 
|| 1x0 1| f[i/2 < fy/Co so this method also provides a proof of Kato's theorem [32]. 
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To simplify notations, let us denote by 5*(t) = mm{<5 (t), Co<5o(r)}. On [0,T] x T 3 , the inequality (62) 
reads : 



Vol(' 



— y 

(J3) 



e s *W\u(t,k)\ 2 < C\\u \\ 2 L2 



fcex 3 



and (9) gives a sequence of dimensionless coefficients ck > such that : 

= (27r)- 2 p( I pj e- 5 * K e 5 *l fc l|u(t,fc)| 2 <CVol(T 3 )r^e- 5 ^^(0)-c^ with y CK <l. 

(64) 

On [0, T]xE 3 , the inequality (62) does not directly provide a pointwise estimate of E*(t, K) so one considers 
instead : 

1 f 2K 

E*(K,t) = — / E*(t,k)dk 
K Jk 

The spectrum is given by (4). One gets : 

• r2K e -6'\t\ 



E*(K,t) = (2tt) 



-3 P 



One checks easily, using the properties of the cut-off function tp (uniform in the limit ij){r) — 1 5 r =i), that 



/■2JT -5*\H\ . /-lei/* x 

/ ^— y>(^ (IKe-^ x / 4;(r)dr < 2c{K)e—2 5 K 
Jk k \ k J J\e.\/2K 



\Z\/K 



with ^2 ne z c (^ n ^) — 1- Inequality (62) gives : 



3 



S*(A",t) < 2(2vr)- 3 e-2 



c(K) 
IT' 



< CK-*-e-* s * K E{G) -c(K). 



(65) 



This concludes the proof of (55) and of Theorem 12. 



Bibliographical note on the analytic estimate (62). For physical reasons, one was interested only in 
L 2 norms, therefore a pointwise majoration in Fourier variables was sufficient. One should notice however 
the following idea suggested by [49]. Instead of using |£| < |£ — r/| + M, one can rely on the exact formula : 



Va E 



(uv) = {e ia ^u) (e ta ^v) . 



One could therefore replace (57) by a family of ODEs 

Mt)= [ Ve iA «M Q ' v n(r) 2 dr with A Q (t) = V^i - \0 a (t) 
Jo LP 



that would provide uniform bounds on e^ Ucx ^ u(t) for t £ [0, T*\ and a £ B = {(3 £ R 3 ; |/3| < C}. One 
can then deduce the local-in-time L p -analytic estimate of u because an elementary exercise in Littlewood- 
Paley theory states that : 



< C sup 

L p a&B 



ut ia-V 



LP 



However, due to the lack of a uniform a-priori bound of ||w(t)|| iP when p / 2, one can only propagate the 
local estimate along the time-line by stating that 



u is smooth, namely u £ L°°([0, T] x 17) e^-™ n ^ vt ^ A u{t) £ L p 



(66) 



but unlike (62), the corresponding inequality involves both uq and sup[ 0j T] 



\Lp(n) 



to compute i* and 



the L p -analytic norm. Statements similar to (66) can be found e.g. in [47], [49], [48] or [59]. However, as 
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these authors make the wise choice not to state the corresponding estimate, the status of (62) was unclear. 
In particular, estimates hidden in (66) cannot usually be translated in time with uniform constants. For 
the convenience of the reader I decided to provide an elementary proof of (62). The method was greatly 
inspired by [18], where a statement similar to (61) was proved for the L 2 ([0,T] x M 3 )-norm of e^^utt). 

Some subtle mathematical problems are closely related to Theorem 12. Analyticity in time and Gevrey 
classes have been studied in [28]. The difficult question of analyticity for general domains seems to have 
been studied only once, in [55], and the question of external forces is dealt with in [34]. Readers interested 
in analyticity questions for 2D-turbulence should refer e.g. to [4]. 

5 Necessary conditions satisfied by turbulent flows 

In this section, one investigates the necessary conditions satisfied by K41-turbulent flows. One proves 
the relation between Integral and Taylor Scale Reynolds numbers (Theorem 14). On T 3 , this relation is 
shown to characterize a time-scale on which free turbulence can be observed (Theorem 15). Finaly one 
investigates whether smooth solutions can be K41-turbulent (Theorems 16 and 17). 

5.1 Reynolds numbers 

(k \ 4 / 3 

The Reynolds number 3ft = (7^) is the so-called Integral Scale Reynolds number [31, §7.18 p. 107]. 
In experiments, one uses mostly the Taylor Scale Reynolds number Ur ™ XT with U Iuls (t) = y/E{t) and the 
Taylor scale defined by t-oVtt = ^J"^! 1,2 ; in other words : 



(*) " IK*)II 2 

TT \ „l/2 IL./MMI 

L 2 



£/rmsA T P 1/2 \\u(t)\\ 2 r 



v\\Vu{t)\\ L2 

Usually R\ is "tailored" to the needs of each experiment to get a time-independent number. 
Definition If u is a Leray solution of (1), the Taylor-Scale Reynolds number of u on [Tq, T±] is : 



Observations [31, 7.17 p. 107] suggest that R\ ~ 3ft 1//2 . This is indeed a rigorous fact. 

Theorem 14 Assume that u is a Leray solution of (1) and a K41-function. 
1. For fl = M 3 , assume also that uq £ L 1 n L 2 . Then : 

3ft~ 97 / 4 jft 163ft 97 / 4 

216(1 -K- 1 )(l-K- 1 /2)2 ~rJ- 27(l-R- 1 )(l-K- 1 /2)2- ( > 



In particular : 



a 3 z;3ft a^Vol^lTcTi]) 1 ^ a 3/2 u 



According to (32), the symbol ~ means that both inequalities hold with constants depending on 3ft 
and 3ft 7 but that the constants have a purely numerical limit in the asymptotic (31). 



2. For T», let US sssume t M t f ^ x)il = M d U (39) fa Sat fa&d. Tien similar pities hold. 

(1-C(n_)) 2 3ft 2048 



't 3 

The numerical values are : 



+ o(^ 1 )<^< — +o(M-, 7 ) (70) 



864 v ' ,; - R\ ~ 3375 



where o(3ft;7) — > in the asymptotic (31). 
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Note that the right-hand side of (70) requires only a lower bound on K_/K c , thus it holds regardless of (39). 
Proof. Using the definition of 3£ and (35), (36a) for Kd and K c , one gets : 

\K-J \K C J o?ve A 

thus Theorem 6 implies 3? ~ R\. Proposition 11 provides Vol(it; [To,Ti]) ~ (K-)~ 3 hence 

3? ~ (Vol(u; [T^T^K^ = ^ Vol(«; [Tq,^]) 4 / 9 . 

Using 3? ~ R\ provides 

E 2 \ 3/4 E 3 / 2 

a 2 Vol(u; [T ,r!]) 4 /9j = a 3/2 Vol ( M; pzj^jji/s" 

Substitution of this formula in the previous expression of 3? yields the last formula in (69). In the case 
of SI = T 3 , the proof is similar. ■ 

5.2 Time-scale on which free turbulence can be observed 

Definition Given a Leray solution u of (1) on S7 = T 3 with / puo = 0, let us define the following time 
scale, that one could call for example the "transfer time from uq to ui = rotu on [To,Ti]" : 

f Tl 2 

n 3„2 / Mt)\\hdt 

T(u ;co) = x JT « (71) 

3 

According to (19b), one has T(uq,uj) < q ^ L - The following statement proves that 3?T(iio,w) is the precise 
time-scale on which turbulence can be observed : for a two short observation time one will not see K41- 
properties and for a too long observation time the time-average will describe the fluid as being mostly at 
rest. This time-scale appears also naturally in the computation of e for a smooth flow (see (77) below). 

Theorem 15 If u is a K41-function on [To, T{\ x T 3 and a Leray solution of (1) with / puo(x)dx = 0, the 
following inequality holds for numerical factors Cj(3t, 7) — > 1 in the asymptotic (31) : 

3375 rf3? ^»rr„ -,a<t t < 128C ^,7) 1 Vq1(t 3 ) 4 / 3 

—^(3?, 7 ) 3? T(u , u) < T X - T < ^T^u) * ~ 2 (72) 

Moreover, if (39) is satisfied, then : 

*™<Wt,y) < (i-)'£z«< < (73, 



2048 dV ' " ~ \E J 3?T(u ;w) ~ (1-C(n_)) 2 ' 

In particular, (72) gives : 



3i 512C 5 (3?, 7 ) .., w Vol(T 3 ) 2 / 3 

^- < „„„ 2 with 3i c = — -• 74) 

3t c 3375^^ uT(uq;oj) 



Proof. Poincare's inequality reads : 



u(t) — I u(t, x)dx 
n 



Vol(T 3 )V3 
< \J Vtx t) L3(T3) 



29 



On T 3 , the total impulsion is preserved so (43) holds for any t > T> Combined with (19) one gets : 

K^lli" + 2 - (^oiWTi) 2 £ IWH* dt ' ^ W< T ^ 

and Gronwall's lemma provides (see [38] for a numerical confirmation) : 

IKt)|| L2(T3) < ||n(T )|| i2(T3) exp (- ( Vq1( ^ )1/3 ) ~ TqU ■ 



In turn, this gives : 



2vr \ 2 ,„ „ A , _ N 1-e 



E[K) - Lu<T>l > Vol(T3)i/3 ) ^- T °)l (|)( >> 



Combining the definition of e and R\ gives : 



Tl „ / m,2 , (Ti -T )T 2 SR 



Theorem 14 on T 3 ensures » < (§§ + o(l))K$, which now reads 

jC im<)& dt ~ (ii + 0(1> ) o^s x (Tl " To)$2 ((wljis)'' 1 ^ " ^ 



and <3? 2 (s) < min{l; 1/s 2 } provides both upper and lower estimates on T\ — To in (72). Note that one does 
not even require (39), for this assumption is only needed to get a more precise upper-bound of T\ —Tq. 

Conversely, if (39) holds, then »>((!- C(n_)) 2 /864 + o(l))i? 2 and (75) boils down to (73). 



Compatibility with the time scale on which smoothness is guaranteed. Mathematicians can 
guarantee the smoothness of the solution of (1) on at least [To, To + 0] with e.g. 

II w o|Il2 

Is such an interval long enough for the observation of free turbulence ? According to Theorem 15, the answer 
is yes provided 6 > $t.T(uo,oj). Indeed, if Tl — To < O < CT(uq,uj) then the only possible Reynolds 
number on [To,Ti] is 3? < C and the turbulent asymptotic (31) cannot be achieved. 

Inequality 6 > ~RT{uq,uj) is equivalent to 

/•Ti 

ii n4 \ Str\ ii n4 / ii ii 2 

P v IPo|| L 2 <, « ^ ||wo||x,2 / IMIl 2 ) 

J To 

which, combined with Poincare's inequality and p = Vol(Sl) -1 , implies : 

— (Ti-To^^a 3 ^ /" Tl ||a;||i 2 < ^Vol^) 1 ^. (76 ) 
P v J T 

For large initial data, this makes it impossible for the initial energy to dissipate almost completely on [To, Ti] 
because the left-hand side would then be equivalent to a 3 $l(pv)~ 1 Eo » ^Vol(Q) 1 / 3 . 

To put it simply, the time-scale on which regularity is guaranteed is too short to observe a fully 
developed free turbulence and thus it might not have a deep physical meaning. One should however refrain 
from jumping to the conclusion that K41-turbulence is an obstruction to smoothness. If one keeps (53) in 
mind (which is accepted in every physics textbook), it indicates instead that (1) can develop a very specific 
dynamic in Fourier space and that the techniques used to prove local smoothness have failed to capture it. 
We will see in §6.8, that K41-turbulence is troublingly close to the best known local smoothing effect. 
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5.3 Two necessary conditions satisfied by smooth turbulent flows 

The goal of this last section is to investigate the necessary conditions that occur when a smooth solution u 
of (1) happens to be a K41-function. Two Theorems can be stated. 

5.3.1 Necessity of intermittency 

In section §1.4, temporal intermittency was defined as a substantial deviation between e(t) and e. If u 
is smooth on [To,Ti], one can compute e with (24) : 

_ _ Eq — Ei 
£ ~ T\ - T ' 

The conservation of energy reads : 

- .. El 

'T 

hence 



E Q -E 1 =pv [ 1 \\oj{t)\\ 2 L 2dt= ^§-T{u -oo) 

JTo a V 



e-=#~- (77) 



Theorem 16 If u is a smooth solution of (1) i.e. u £ L°°([To,Ti] x Q) with = R 3 or T 3 and a 
K41-function on [Tq,Ti], the following condition must be satished : 



^ ifn = 



2E n .L E n ~ V ' U \\ T-y-T Q | lSyTs ifn = T 3 

with C(§R, 7) — >■ 1 in the asymptotic (31). The numerical constant is bounded from below by 1.299 on 
and by 1.283 on T 3 . 

Proof. Theorem 14 provides in the asymptotic (31) : 

^_ E 2 ft < E 12 +o(l) ifft = R 3 , 

« 3 Z^ X i? 2 ^ ^Sfl X I 2048 + o(1) iffi = T 3 



Estimate (25) reads : 



^ 2 < (E + E 1 



+ y |e(t)-e|d^ 



I To 

and (78) follows immediately from the comparison with (77). ■ 

One can wonder whether (78) actually provides a lower bound on intermittency. For SI = T 3 , the answer 
is subtle but positive. According to Theorem 15 and a careful track of constants, the right-hand side of (78) 
cannot asymptotically exceed 1. On the other hand, the energy estimate gives E< 2Eq 1 £ [1/2; 1]- Therefore, 
if one assumes that : 

• most of the initial energy has been dissipated, i.e. E\ <C En, then E °^ 1 — \ 

• and that the observation time is the minimal, i.e. T\ is chosen short enough such that the right-hand 
side of (78) belongs to ]\, 1] 
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then the left-hand side can be bootstrapped in the right-hand side and inequality (78) indeed provides a 
lower bound of intermittency : 

I 1 \e{t)-e\dt > ^E 

with C G]0, 1]. Given smooth data uq, it is not clear whether one can find T\ such that both conditions 
are simultaneously satisfied (the problem is to prove that E{Tq + $lT(uo; oj)) < cE(Tq) with a sufficiently 
small numerical constant c < 1). However, physical intuition suggests that this is the case since the best 
time-scale to observe free turbulence is to wait till most of the initial energy has been dissipated but not 
any longer. 

5.3.2 Compatibility of the K ~ 5 / 3 law with the analytic smoothing effect 

The question dealt with in this section is the following : are the finer scales of turbulent vortex structures 
limited by the analyticity radius of the solution ? In other words, if 8 denotes the analyticity radius of a 
K41-turbulent solution u, what is the possible range of 5K + ? In the regime 5K + < C the finer scale 
of turbulent structures is K7 1 and is limited from below by C~ 1 5, which drastically limits the possible 
Reynolds numbers : 

/ C \ 4/3 ( C L\ 4/3 

3ft < I ) and in particular < I ) on T 3 . 

\oa_ / \2ir o J 

Conversely, the regime 5K+ S> 1 means that turbulent structures exist at much finer scales than the 
analyticity radius and the asymptotic (31) remains possible (see also §6.7.2 below). 

The following result denotes the compatibility at K = K + between K41-property (33) and the high- 
frequency bound on the spectrum given by Theorem 12. There are two cases depending on whether the 
initial data is supposed to have no additional smoothness to being L°° n H 1 ^) or if on the contrary, one 
considers a flow already "well prepared" by the analytic smoothing effect. 

Theorem 17 Let u be a smooth solution of (1) on [T ,Ti] x with Q = R 3 or T 3 and a K41-function 
on [Tq,Ti]. There exists a numerical constant C > such that the following conditions are satisfied. 

1. Unprepared data on T 3 - If u(T ) G H 1 ^ 3 ) with / pu (x)dx = and if u £ L°°([T ,Ti] x Q), 

J ■■■■ 

then : 

j < C ^ [eT**+ + 6 ° K+ ^L_\ (79) 
K 3 7/2 E(T ) ~ CK f + 1 + (5 K + ) 3 Tx - T / { ^ } 

C r- C v 2 
2 2sup [TojTl] ||a;|| L 2 

by Theorem 12 on [To + r, T±\. 



provided (39) holds and where So = —V^t = — ~2~ is the analyticity radius guaranteed 



2. Well prepared data on T 3 — If moreover one assumes that the initial data uo £ i? x (T 3 ) was given at 
t = and that the K41-property holds on [To, T±] with 



v 3 



To > (80) 

S U P[0,T!] IM| L 2 

Co 1/2 

then Theorem 12 guarantees an analyticity radius of at least 5\ = ^~ on [To, T\\. In this 

2 sup [0jTl] || Co- 1| 

case, (79) can be improved to read : 

{ ' K ^ log (T^)) +log (T) +21ogR ( 81 » 
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3. Case of R 3 — Similar results are also valid for = M 3 except that the assumptions J" T3 puo(x)dx = 
and (39) must be dropped and replaced by uq G L 1 n L 2 (R 3 ). Estimate (79) reads : 

l -^' 2 " <C^\ E -^ + ^^^) (82) 



(83) 



^37/2 £( To ) " \ 1 + (5 K+)3 Ti - T 

and (81) becomes instead : 



Remarks 



1. According to the comments at the end of §5.2, the "well-preparedness" assumption (80) is likely to be 
satisfied if the initial data is large enough. 

2. In the unprepared case, one almost gets the asymptotic 5qK + < (E /E) 3? on T 3 and (Eq/E) 1 / 2 3ft 1 / 4 
on M 3 . In the well prepared case, it improves rigorously to : 

* lK+ <C< + lo g(£o /£) + {f gR ° nn = t ' 



2 



log 5ft on Q 



It will be shown in §6.7.2 that there is an experimental hint towards 5K + > 1 and that 5K + ~ 1 must 
hold when physicists claim to observe fully developed turbulence. 

3. One can compute 5iK + {i = 0, 1) using only norms of the vorticity : 

^5/4^1/4 

6lK+ ~ «3/4 supi ||a;(t)||i 2 ' 
Thus, using E 2 ~ q 3 i/eS from (69), the well prepared case can be reformulated in the following form : 



CI/ 5/4^i/4 \ _ 9 [K 3 on n = T 3 

TT 



o?ve ■ exp — =- < C£ 2 x <( "~ "" (84) 



Proof. Let us first investigate the case of T 3 . Comparison of Theorem 12 with (33) for K = K + provides : 

3T 37/4 ae 2 / 3 i^ 5/3 < E*{K+) < CK + L 2 (e' 5 ^ K +) E(T ) 
where S(t) is defined in Theorem 12. Let us write it down as : 

3r 3 ^ 4 a£ 2 / 3 K; 2/3 < C(K + L) 2 (e- 5 ^ K +) E(T ). 
According to definition (36b) and Theorem 6, one has : 



ae 



■2/3^-2/3 = (KA 2 "e > 1 " C ^ JL 
' K + ) ~ 12K 3 7/4 Rl/2 



Similarly, (K + L) 2 = (K_L) 2 K 3 / 2 < (27r) 2 3? 3 / 2 . A direct computation of the time-average of e S( $ K + reads : 

e-^ + )=e-^ + + _^_$(^ +) 



± ._ JC Co 

II 1 1 4 

su P[T ,Ti] IMIl 2 



with r = — ^— and <5q = — \[vt and the function 



«(.) = C§ (' 2 "-" + ^' - - e- 
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1 / C 2 s \ ( C 2 s \ 

that satisfy - I ) < VK S ) < 3 ( °J for any s G R+. This gives (79). In the case of "well 

prepared" data, Theorem 12 gives S(t) = 5o for any t G [To,Ti] thus (e~ s ^ K ) = e~ s ° K+ and 

l-C(n-) _g_ 2 ^ + 

5R37/2 E(0) - UJt e 

from which (81) follows immediately. 

On M 3 and provided mo 6 L 1 , a comparison of Theorem 12 and (33) implies instead : 



which explains the different game of powers of the Reynolds number -ft. 



6 Final remarks and some open problems 

Let us end this article with some comments on the questions at stake and especially the compatibility 
between the spectral theory presented above and the experimentally accessible quantities called structure 
functions. I will conclude with a striking numerical fact that suggests that turbulent flows might actually 
saturate the current estimation of analytic regularity. 



6.1 Finding examples of turbulent flows 

The burning question is to find examples of K41-turbulent flows. According to our definition, this 
problem can be split in two. The first step is to find K41-functions that are solutions of (1). The second 
one is to find, among those functions, some that satisfy the asymptotic (31). 

Proving that a given function is K41 requires to compute a lower bound on Vol(u; [To,Ti]) and then 
check that a substantial amount of enstrophy is contained in the scales below this bound. It can at least be 
tested numerically as this property is stable in Leray space. Satisfying (31) is more subtle and is likely to 
require additional assumptions. The most obvious one will be a form of isotropy (with a proper definition, 
still to be found) because the energy spectrum was defined with isotropic spectral cutoffs. The question 
of an anisotropic theory is also widely open and could be of interest in concrete situations. For example, 
atmospheric turbulence in the jet stream (rapid air flows at high altitude, well known to jet pilots) is 
anisotropic because the ratio height/width is a tiny parameter. 

For bounded domains, one could for example investigate solutions whose vorticity satisfies : 

oo(t, x) 



n |w(t,s)| dX 



< c - 



The cheapest conjecture is that for Co small enough (which means that the vortex lines are somehow dis- 
tributed isotropically), the asymptotic (31) should hold. 

This is also where probabilities might prove handy. The corresponding conjecture is that (31) holds in 
an average sense over a statistical ensemble of solutions, for some ergodic probability measure. 

However, these conjectures should be balanced by taking into account substantial fluctuations of the 
dissipation rate e both in the temporal and in the spatial domain, i.e. intermittency (see §6.2). 



6.2 Geometrical structures of turbulence, localization and intermittency 

Understanding the geometrical structures involved in turbulence is a major challenge and the core of 
modern research on turbulence (see e.g. [9], [22], [63] and the references therein). The main idea is that 
substantial fluctuations of the dissipation rate exist both at large scales (caused by the mechanisms of 
agitation) and at small scales (caused by the stretching of vortex lines). In consequence, a refined theory of 
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turbulence cannot rely solely on the average value e. In this article the focus was on "large scale" turbulence, 
i.e. in the case where the fluctuations of e are small compared to e. 

The properties of Vol(u; [7o,Ti]) and T(uo;uj) in regard to the characteristic scales of the geometric 
structures in turbulent flows call for a closer look. For example, a starting point for studying intermittency 
is the following definition. 

Definition Given Cl'(t) C Q a smooth family of smooth subsets of 0, a function u G £(fi) is said to be a 
local K41 -function on Q' x [Tq,T\] if 

u'(t,x) = u(t,x)x(t,x) 
is a K41-function on x [To,Ti] where x is a smooth cutoff function such that 

X (t,x) = 1 if x € n'(t) and Vol( X ; [T ,Ti]) ~ Vol(l n /; [T ,Ti]). 

The first question of localization and intermittency is to study which sub-domains of Qq x [To,Ti] are 
admissible if u is a local K41-function on f2o x [Tb,Ti] and to determine the corresponding parameters 
K±(SV), and 7(f)') • Using (37d), an obvious restriction is Vol(u') < Vol(x) < Vol(u). 

Related to the volume function is the following inequality, which (with Leray's inequality) is one of the 
only estimates that has a sub-linear growth in time : 

\\u{t)\\ L i {m3) < \\u Q \\ L i + \\u \\ 2 L 2 . 

Improving this inequality for a given flow on R 3 or finding examples of saturations on some time intervals 
would definitely be interesting. 

6.3 Modelisation versus PDEs 

When dealing with the spectral description of turbulence, two natural questions occur. Which solutions 
of (1) satisfy the spectral asymptotic (31) ? And in that case, what are their qualitative properties ? 

Historically, Kolmogorov addressed both questions simultaneously by proposing a model based on strong 
probabilistic assumptions that cannot be checked directly, namely by assuming that "the difference u(t, x) — 
u(t',x') is a probabilistic variable whose law does not depend on t,t',x,x' and is not affected either by 
rotations of the coordinate system" [40] . Obviously the consequences of his predictions have been extensively 
studied and most of them where found valid, but sometimes with not as good a precision as one could hope 
for : structure functions S p (£) (see §6.7) are predicted to scale as £ p ^ 3 but for p > 4 one observes that 
the exponent p/3 must be corrected by some small negative term. These mishaps are known under the 
generic name of "intermittency" and various models have been sought to explain them, including a second 
probabilistic model by Kolmogorov [43]. 

The main contribution of this article is to show that one can study the qualitative properties of turbulent 
Bows independently from the a-priori models of the structure of such a Bow and that it can be done using 
deterministic tools. This path should allow colleagues to concentrate on the core mathematical problems. 

6.4 Spectral problems for general domains 

One cannot deny that the tradition of probabilistic models in turbulence is a convenient way to avoid 
dealing with the spectrum of the Stokes operator on domains. . . However, the lack of precision on the 
true energy spectrum (the function denoted by E\ defined on a {A 1 / 2 ) where A is the Stokes operator) 
considerably darkens the foundations of some experimental protocols that focus on the spectrum E* defined 
with Fourier coefficients and assumed without proof to be the shell averages of . 

More precisely, as the spectral theory of the Stokes operator is not known, the naive protocol is to collect 
data in some subregion, then compute Fourier coefficients as if the flow was periodic, using an FFT-type 
algorithm with ad-hoc anti-aliasing techniques, and then finally compute the energy spectrum with the 
formula (13) valid for T 3 . The question is : does one really look at the spectrum of u ? 



35 



For example, sorting the spectrum of the Stokes operator on ]^[(]R/LjZ) and checking (14) raises non- 
trivial problems of rational approximation as soon as Li 7^ Lj. However, (14) played a crucial role in the 
proof of (13) on T 3 . Worse, if one assumes generically simple eigenvalues, Weyl's asymptotic for the Dirichlet 
problem in dimension n gives K? ~ C(n; \Q\)j 2 / n which suggests 

K K < C'{n-M) 



with no general lower bound if the domain is a small perturbation of one with multiple eigenvalues. 

Of course, other less naive experimental protocols are used but they also have their share of mathematical 
problems (see §6.7.2 and §6.7.3 on structure functions). 

The author conjectures that the properties of E^{K, t) : a{A 1 / 2 ) x R + — > R + reflect a general arithmetic 
of the specific nonlinearity ¥((u ■ V)u) and thus could be truly universal. However, discrepancies between 
the various possible spectra of the Stokes operator might explain why E* defined by (13) is not always the 
appropriate experimental quantity. 

This might be a first key explaining why there seem to be "many" theories of turbulence, depending on 
what kind of flow one deals with ; it is likely that the spectral theory for the Stokes operator in an outer 
domain has little in common with the spectral theory in a wind tunnel or that of a swimming pool. . . 



6.5 Universality of the Kolmogorov constant 

Conversely, to give some credit to the protocols based on computing the sum of Fourier coefficients on 
spherical shells of frequencies, one could investigate the properties of sub-domains of turbulent flows and 
show that "away from anything" the rule of thumb is that of M 3 or T 3 . 

For example, experimental evidence suggests that the Kolmogorov constant a in (33) is universal and 
that, with reasonable precision, it does not even depend on the flow or the shape of the domain (provided 
turbulence is homogeneous and isotropic). Published values are in the range a G [0.45,2.4] with a common 
agreement [62] around a = 0.5. See also [67], [36] and the references therein. The proof of the universality 
of a small range for a would be very instructive 3 . 



6.6 Scale-by-scale balance of energy and energy cascade 

The following identity describes the scale-by-scale energy budget and follows directly from an energy 
estimate of (1). Let us denote by Sk = X * ne l° w -P ass filter and S*^ its adjoint on L 2 ($7) : 

s(f *"**>*Hs ('■*«*) 

= -2pu \\V(S K u)\\ 2 L2 - 2p (S* K S K ((u • V)u)\u) L 2 
+ 2p ([S* K S K , uA]u - [S* K S K , V]p\u) L2 
On M 3 or T 3 , all the commutators vanish and one gets : 

j t ^ K E*(k,t)dk^j =-2pv\\V{S K u)\\ 2 L2 -2p{S* K S K {u®u)\Vu) L2 . (85) 

This identity is the base of the famous "energy cascade" interpretation. The left-hand side denotes the rate 
of change of the energy contained at low frequencies, i.e. "at scales larger than K -1 " . It is balanced by 
the energy dissipation at such scales (first term on the right-hand side) plus the so called "energy flux to 
smaller scales" (second term). If this term is negative, energy is indeed transferred away from [0, K] hence 
goes towards lower scales. There is no "energy injection" term as in our case the external force is zero. 

3 As a "per unit of mass" theory, the transformation (u,p,p) M> (u,fxp,fj,p) was disregarded before. However, to the best of 
my ignorance, the following experiment has never been done : compare the Kolmogorov constant a for flows of fluids having 
similar viscosities but a fixed mass ratio between molecules, when the the same number of molecules and similar velocity fields 
are involved. For example, a water flow and a flow of heavy water where hydrogen atoms are replaced with its heavier isotope 
deuterium. Would the Kolmogorov constant be identical ? 
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Definition A solution of (1) on SI = M 3 or T 3 has the "energy cascade" property at time t > if there 



exists K*_ < K\ such that : 



VK G [K*_,K* + ], V / S* K S K ( Ui ® «j) • $Uj(t, x)dx > 0. (86) 

In the limit of non-smooth frequency cutoff, S* k Sk — > J 7 ^ 1 o 1\£\<k o T . 

There is yet no mathematical evidence that this property holds for a general class of solutions of (1), 
not even for the few explicit solutions known. A rigorous connection between this definition of the energy 
cascade and K41-turbulence is an open problem, in particular the relation between K± and K±. 

Let us however mention this interesting property that was pointed out to the author by [58] . 

Theorem 18 Assume that u is a Leray solution on R 3 and that for all t 6 [To,Ti], one has a "reverse" 
cascade : 

VK>0, V / S* K S K (ui <g> uj) ■ diUj(t, x)dx < 0. (87) 

Then the energy equality \\u(T )\\ 2 L2 = ||«(Ti)||^ 2 + 2v I \\Vu{t)\\ 2 h2 dt holds. 

JTo 

Proof. For all K > 0, assumption (87) and the identity (85) provide : 

\\S K u(T )\\ 2 L2 <\\S K u(T )\\ 2 L2 -2 [ \s* K S K {u®u)\Vu) L2 = \\S K u{T l )\\ 2 L2 + 2v f ' \\V(S K u)\\ 2 L2 . 

JTo JTo 

Letting K — > oo, one gets : 

||^T )||| 2 < \\u{T l )\\ 2 L2 +2v \\Vu\\ 2 L2 < \\u(T )\\ 2 L2 

JTo 



pTi 
I To 

the right-hand side being the classical Leray inequality. 



6.7 Structure functions S p (£) 

As the energy spectrum is not always available, other universal properties have been looked for as 
substitutes to (31) and are commonly accepted as experimental evidence of turbulence. However, the 
connection between these properties and the spectral definition of K41-turbulence is not as rigorous as one 
could hope for and constitutes an immediate source of interesting mathematical problems. 

6.7.1 The second-order structure function 
Definition Let us introduce the correlation function : 



T(t, y) = p u(t, x + y)u(t, x) dx = E(t) - S 2 (t, y) 
JKi 



(88) 



with 



S2{t,y) = ^l \u(t,x + y)-u(t,x)\ 2 dx. (89) 

2 JR3 

The following quantity is called the second-order structure function of u on [Tq, T±] : 

V£>0 S 2 {£) = [ S 2 (£9)d9. 
Recall that the S 2 (y) denotes the time average (23) of S 2 (t,y) on [Tq,Ti]. 



(90) 
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The following result is sometimes called the Wiener-Khinchin formula. 
Theorem 19 For any function u G L 1 ([T , Ti]; L 2 (R 3 )), one has : 

SM _ B r ^m^ K)dK -sp-, SHK)iK , m 



4tt J IK v ' J V ^ 

It is usually suggested in physics textbooks that "the energy spectrum is the Fourier transform of the 
correlation function" but the formula is systematically left unstated. As far as proof is concerned, it is 
usually claimed to be a consequence of various probabilistic assumptions. One can prove instead that it is a 
perfectly determinist statement that relies on the following property : the Fourier transform in R 3 commutes 
with the process of replacing a given function by the radial one whose values are the averages of the initial 
function on each sphere. 



Proposition 20 For a function f G <S(R ), let us define : 

F(r) = [ f(r9)d0 and G(A) = f f{X9)d6 (92) 

J§2 J§ 2 



and S(a) = a sin a. Then, one has 



The inversion formula reads : 



[•CO 

\ 2 G(X) = 4tt / F(r) S(r\)dr. (93a) 
Jo 



r 2 F(r) = —j / G(X) S{r\)d\. (93b) 
2?r Jo 

Moreover, if f(x) = ^F(\x\) is radial, then f is also a radial function thus /(£) = |^G(|£|) can be computed 
with (93a)-(93b). 

The radial statement is classical (see [64, Chap. 4, Theorem 3.3]). The rest is implicit to the stability under 
Fourier transform of the orthogonal decomposition of L 2 (M 3 ) into spherical harmonics [64, Chap. 4, Lemma 
2.18]. Indeed, joining I? = S)j and T : S)j — ^ S)j implies that T commutes with the orthogonal projection 
on each fjj. This abstract argument allows one to claim that ^F is defined for any / G L 2 (M 3 ) as its 
orthogonal projection on S)o and that the statement still holds in this case. 

Proof. Let us first establish the statement for a radial function /(|x|) on R 3 . Its Fourier transform is : 

\tfH0 = J™ s(r\t\;-^f(r)dr with S(p,0') = p 2 'dB. 

The function S : R+ x S 2 — > C is invariant by rotations of the second variable : 

S(p,9) = S(p, ei ). 

As p^ 1 S(p, e±) has the same value and derivative at the origin as Ansinp, one has S{p,e\) = 47rpsinp (one 
could also directly compute the integral in polar coordinates) and 



/•OD 

7(0=4*/ f(r) S(r\£\)dr. 
Jo 



This proves the theorem in the case of a radial function. 

Let us now turn to the general case. One defines a function /o by : 

f(x) = ±F(\x\) + Mx) 
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which ensures that : 



Vr > 0, / f o (r0)dO = 0. 

J§ 2 

Applying the theorem for the radial part gives : 

G ( A ) = S r F(r) S{r\)dr + f f (\9)d9. 
A Jo Js 2 

The last term boils down easily using Fubini's theorem : 

/ fo(\0)d0 = ^ // e-^fo (\ ) dydO 
Jn 2 A J JRtxs 2 VA/ 



" A3 //(_ 



~ lp9 - e 'h {{e) P 2 dpd6de' 



A 3 IIr+xS 2 P ju ^A 
= 0. 



smp fo^e)p>d P de 



The second to last identity is the previous computation of S(p,9) = 47rpshi/9 and the last one is the fact 
that the sphere averages of /o vanish. ■ 

Proof of Theorem 19. The first step is : 

S 2 (£) =4ttE- [ Y{£9)d9 
Js 2 

which results immediately from the definitions. Next, one observes that : 

p\u(t,0\ 2 = [ e- i ^F(t,y)dy (94) 

which express the duality between convolution and regular multiplication. Taking the time average on [To, Ti], 
then sphere averages, and finally substituting the definition (5), one gets : 

jj(K9)d9= i ^&{K). 
The conclusion now follows directly from Proposition 20. ■ 

6.7.2 Range of validity of the "2/3 law" and radius of analyticity of u 

Experimental evidence [31, p. 57-61] suggests that : 

We Mo], S 2 (£) ^ (3(e£) 2/3 (95) 

where r\ ~ K^ 1 is the dissipation scale and £q — KZ 1 is the size of large eddies. Usually, physics courses 
state that this so called 2/3 law is equivalent to the K~ 5 ^ 3 decay of the energy spectrum. This "equivalence" 
calls for a closer look. 

The poor man's argument is the following. Let us consider an ideal case where u would be a radial power 
function whose energy spectrum is exactly K~ 5 ^ 3 (this means = |£| -11 ^ 6 even though this function 

is not a solution of (1) nor even a square integrable one). Then, for any £ > and 9 £ S 2 , (94) gives : 

f(£9) = 3 T(9). 
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Figure 5: Numerical investigation of the domain of validity of (95). 

Top four - 52 CO is computed by Theorem 19 using an idealized spectrum E$(K) = K~ 5 / 3 l\i m (K) for different 
values of R. Plot in Log-Log scale on [fl^.lO]. Also represented is (R£) 2 S 2 {R~ 1 ) for £ < 10/R and £ 2 / 3 . 
One observes a close fit of £2 CO and £ 2 I 3 on [10/R, 1]. The graph was obtained by formal integration with 
Mathematics©. 

Bottom left - S2CO computed for a "real" energy spectrum 



E\K) 



(K 2 K<1 
K- 5 / 3 Ke[l,R] 

R - 5 /3 e -S(K-R) R>R 



with R = 10 3 and a radius of analyticity 5 G {10 -2 , 10 -3 , 10 -4 } (this time with numerical integration). The 
range of validity of (95) is maximal for R5 = 1 but drops drastically when R5 -C 1. 

Bottom right - Spectral precision ^log5 2 C0 - I for an idealized spectrum and R G {10 2 , 10 3 , 10 4 }. When 
the precision is in the gray band, (95) is satisfied with a relative error of less than 10%. Also represented is 
the precision for the previous "real" spectrum with R = 10 3 and 5 = 10~ 3 . For £ > 10, the function S2CO is 
oscillatory, which reduces the precision of the numerical integration in the case of the "real" spectrum. 
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However in this case E = +00 and 6*2 (^) are undefined. Even excellent physics textbooks do not further 
justify the "equivalence" between the 2/3 law and the K~ 5 / 3 except by a vague reference to a probabilistic 
version of Theorem 19. . . 

Slightly more careful physicists [31, p. 87] state that the 2/3 law is only an asymptotic property that 
holds provided the limits are taken in the proper order : T\ — To — > 00, then v — > 0, then £ — > 0, and that 
taking the limits in any other order will lead to trouble. From a mathematical point of view, a rigorous 
upper bound for S^^) can be found in [23]. However, the relation with (31) is not established. 

One can check numerically that (31) and the 2/3 law indeed share a strong connection (see Figure 5). 
The conclusion of this computation is that the domain of validity of (95) is [CKZ 1 , KZ 1 } if the analyticity 
radius of u exceeds KZ 1 - However, if the analyticity radius becomes smaller than KZ 1 , then the range of 
validity of (95) shrinks dramatically. 

This fact should be put in perspective with a common experimental fact : physicists observe the range 
of validity of (95) to be often of much smaller amplitude than the inertial range. Figure 5 suggests that the 
analyticity radius S of such a flow is smaller than Kolmogorov's dissipation scale KZ 1 - Conversely, a fully 
developed turbulence for the structure function S2 {£) indicates that 5 > KZ 1 . 

Even though this numerical study is encouraging, it leaves the rigorous connexion between K41-turbulence 
and the structure function S^) on the list of open problems. Moreover, as the numerical observation sug- 
gests that S2CO = 0(£ 2 ) around I — > thus one cannot expect to prove (95) by a finite expansion of S^^). 
Instead one will have to prove directly that there exists a smooth function 70 with 70 (s) <C 1 if s > 1 and 
such that 



£e[CKZ L ;K: 

where 5 denotes the radius of analyticity of a K41-function u. This subtle exercise in harmonic analysis 
should be an excellent warm-up round before tackling the question of finding examples of turbulent flows. 
This also explains why physicists dodge the problem of accessing the "intermediary" regime of (95) by taking 
suitable asymptotics that will push its domain of validity all the way to £ — > 0. 

6.7.3 Higher order structure functions and the "4/5 law" 

Other structure functions play a central role in experimental protocols. The most celebrated is the so 
called "4/5 law at inertial scales" for the third-order function (this function was historically related to the 
measurement mechanisms used to acquire experimental data in real flows) : 

s " w =</L, (w ^ +w) -" (t ' a;)) ' 9)3 ^)> * -\« (%) 

Recall that the brackets (•) denote time average. This fact is claimed as being "rigorously established" in 
most physics textbooks (see [9, chap. 2] or [31, p. 76-86]). It was indeed addressed in the first paper of 
Kolmogorov [40]. However, the rigorous path from (31) to (96) still requires some enlightenment. Worse, 
one can conjecture that the property Ss(£) = C£ + 0{£ 2 ) holds independently from the asymptotic (31) and 
that it is only a consequence of the rapid decay at infinity of the energy spectrum. Let us follow a proof 
step by step and point out the dark spots. 

The starting point is the following identity, called the Karman-Howarth-Monin relation (let us recall 
that T has been defined by (88)) : 



d t T - 2vA y T = a with a(t,y) = ^ Tr {\u <g> u)(t, x) ■ [Vu(t, x + y) + Vu(t, x - y)]} dx. 

2 J R 3 



(97) 



Even though it is often presented as a consequence of the probabilistic assumptions, it is a perfectly de- 
terministic relation that follows immediately from equation (1), the definition of T and the obvious fact 
that : 



AT(t,y) = Vu(x + y)Vu(x)dx. 
2 Jr3 
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One can rewrite the right-hand side in the following way : 

a(t, y) = ^ dw y \5(t, x, y)\ 2 S(t, x, y) dx^j . (98) 

with 5(t, x, y) = u(t, x + y) — u(t, x). 

Let us now compute the Fourier transform and integrate over a sphere of radius K. One gets : 

d t E\K, t) + 2vK 2 E\K, t) = d K U(K, t). (99) 

One can compute H(K, t) directly from a(t, y) using Proposition 20 : 

U(K,t) = -^[ K [ f e^y d a{t,y)k 2 dkdedy 
1 f sm(K\y\) - K\y\ cos(K\y\) 



7T 2 7r3 \y\ 3 



where a${t,l) = / a(t,W)d6 and / dt = 1. Next, one computes the time averages on [To,Ti]. 

J§2 Jq 7T I 

The Fourier inversion formula in Proposition 20 then reads : 

f f°° sin KP 

(1 + £de)a = — / — — U(K)dK. (100) 

87T J K 

The time average of (99) over [To,Ti] then integrated over [K-,K] reads : 

H(jo - w-) = *> f K ^w* - f Et(K ' T : > "f ( "- T ' ) ^ 

If the spectrum decays sufficiently fast, the right hand-side is almost constant in the range K > K + . Each 
term is also equivalent to e if the solution is smooth. Thus, for K > K + , one gets n(i^) — Il(K-). Usually, 
this claim is the first dark spot justified by some "suitable" asymptotic like letting t — > oo, then v — > with 
fixed e. It seems however to be merely a consequence of the analytic regularity. 

The last step is to use (100) to convert the constancy of lloo = fl(K) = f[(K ) for K > K + . The 
right-hand side can be developed as follows for t — > : 

This singular ODE admits only one bounded solution as t — > 0, namely ao(£) = + 0(£). Substituting 
the definition of a(t,x), one gets : 



/ r,\ id 



2 [[ \5(t,x,£6)\ 2 5(t,x,e9)dxd9 
It is a weak form of (96) that implies by integration of the finite expansion : 



U ^ + 0(£). 
4 



If \5(t,x,£e)\ 2 5(t,x,£6)dxd6 = Cl + 0{£ 2 ) (101) 

for some numerical constant C = ^IIoo. As the remainder terms have been neglected, this identity is not 
rigorously established but one can conjecture that it will hold for any smooth solution of (1) whose spectrum 
decays sufficiently fast. 



42 



Using a probabilistic assumption of homogeneity and isotropy in the region of observation f2, physicists 
claim that : 



Vol(fi) 
96 



(3 + £d t )(5 + £d t ) 



5. 



This identity is the second dark spot because it is not clear how to get a similar formula independently of 
any a-priori model of the flow. Substitution in the previous equation for ctq then gives : 



1 

12 



(l + ede)(3 + £d £ )(5 + 



and this singular ODE admits only one bounded solution as £ — > 0, namely S^(£) = — ^H^i. One can check 
immediately that Iloo has the dimensions of a dissipation and physicists claim indeed (and this is the third 
dark spot) that floo = e. 



6.8 Extremal properties of turbulent flows regarding the optimality of analytic esti- 
mates. 

Let us put an end to this section and this article with a striking observation inspired by a comment of 
Claude Bardos : "Personally I do not believe that the solutions of the incompressible Euler or Navier-Stokes 
equations blow up, but it may well be that there are no other general estimates than the one presently 
found" [1]. 




Left - Plot of xs( K ) in Log-Log scale for 5 = 10" 4 . 
The straigh line is 0.85K~ 5//3 . The upper graph is the 
rigorous bound K~ 1 e~ SIi of the energy spectrum. 
Top - Relative error xs(K)/(0.85K- 5 / 3 ) in Log-Log 
scale for 5 G {KT 1 , 1(T 2 , 1(T 3 , 1(T 4 , 1(T 5 }. For 
5 = 1(T 3 , it shows that Xs(K) matches 0.85K~ 5 / 3 
with less than 25% of relative error on almost 4 decades. 



Figure 6: A striking observation : a logarithmic correction to the rigorous analytic estimates of the energy 
spectrum has a definite K~ 5 / 3 behavior on a large range of K. 



Let us recall our upper bound of the energy spectrum for smooth "old" solutions : 

E*{K) < K- X e- SK E{T). 

This estimate has been proved by a technique that seems to be the cutting edge of quantitative smoothness 
estimates for parabolic equations. It is natural to ask whether flows exist for which this inequality is optimal 
on some large range of K. If this is the case, then the energy spectrum of such flows would exceed any 
substantial correction to this estimate. 

Let us plot therefore the following logarithmic correction (see Figure 6) : 

K -l p -5K 

X5{K) 



\og\2 + K) 
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An extremely troubling fact is that, on a log- log diagram, this corrector shows a definite K~ 5 / 3 behavior ! 
For example when 5 £ [1CU 5 , 10 _1 ], the K~ 5 ^ 3 behavior appears for roughly K G [1, \]. 

This observation is a powerful suggestion that K41-turbulent flows might exist among smooth solutions 
of (1) and that these flows are responsible for the failure of extending local regularity methods. They will 
nonetheless provide examples saturating the classical inequalities of fluid mechanics. 
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